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FOREWORD 


This  work  was  supported  by  the  United  States  Air  Force 
under  Grant  AF0SR-76-3005  and  was  undertaken  at  the 
University  of  Salford  during  the  period  1 June  1976- 
31  May  1979.  The  major  technical  assistance  was  pro- 
vided by  the  research  associate  Dr  D Daintith  but  important 
conceptual  contributions  were  made  by  Dr  A Bradshaw  and 
Professor  J J D'Azzo.  This  final  report  was  expertly 
typed  and  assembled  by  Mrs  D Millward. 
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1.  INTRODUCTION 


There  has  long  been  a crucial  need  for  the  development  of 
techniques  sufficiently  powerful  for  the  routine  computer- 
aided  design  of  on-board  digital  controllers  for  airborne 
systems.  In  order  to  fulfil  this  need,  however,  it  was 

essential  that  the  controllability,  observability,  time- 

« 

optimality,  and  eigenstructure-assignability  properties  of 
multi variable  discrete- time  dynamical  systems  be  clarified 
by  fundamental  system-theoretic  research  since  these  pro- 
perties cannot  be  adequately  characterised  by  classical  z- 
transform  techniques. 

This  report  outlines  both  the  fundamental  system-theoretic 
research  and  the  parallel  development  of  design  techniques 
which  have  lead  to  the  production  of  the  comprehensive 
software  package  EIGENFORTRAC^  . The  use  of  EIGENFORTRAC 
greatly  facilitates  the  design  of  high-performance  multi- 
variable  digital  control  systems  for  a wide  range  of  flight- 
control  applications. 

2.  SYSTEM-THEORETIC  RESEARCH 

2.1  State-Feedback  Regulators  and  Observers 

(2) 

The  results  obtained  by  Porter  (see  Appendix  1)  completely 
characterise  the  entire  range  of  possible  finite  settling- 
time state-feedback  regulators  and  observers  by  specifying 


the  permissible*3*  Smith  canonical  forms*4*  of  the  closed- 
loop  plant  matrices  of  discrete-time  multivariable  systems. 


The  design  of  such  state-feedback  regulators  and  observers 
can  be  readily  effected  by  the  method  of  entire  eigenstructure 
assignment*5*  *6*  and,  in  particular,  by  the  algorithms 
developed  by  Porter  and  D'Azzo*7*  *8*  *9*  (see  Appendices  2, 

3,  and  4) . These  algorithms  readily  yield  the  vectors  which 
are  required  for  the  simultaneous  assignment  of  Jordan 


is  constrained  only  by  the  requirement  that  the  eigenvectors 
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of  the  closed-loop  plant  matrix  must  lie  in  just  one  family 
of  well-defined  subspaces. 
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2 . 3 Dynamic  Compensators 

The  severe  constraints  on  the  closed- loop  eigenstructure 
assignable  by  output  feedback  imply  that  it  is  frequently 
impossible  to  achieve  satisfactory  closed- loop  behaviour 
by  means  of  static  output-feedback  regulators,  and  that  it 
is  consequently  necessary  to  introduce  dynamic  compensators^ 
However,  the  results  obtained  by  Porter  and  Bradshaw ^ 

(see  Appendices  7 and  8)  indicate  that  the  design  of  such 
dynamic  compensators  can  be  effected  by  applying  the  method 
of  entire  eigenstructure  assignment  to  appropriately  augmented^ 
systems.  In  this  way,  the  use  of  observers  can  be  avoided  in 
the  design  of  error-actuated  multivariable  tracking  systems 
even  when  the  special  conditions  previously  established  by 
Bradshaw  and  Porter  (see  Appendix  9)  for  the  existence 

of  such  error- actuated  tracking  systems  are  violated. 

Indeed,  in  view  of  these  fundamental  new  insights  into  the 
structure  of  linear  multivariable  systems,  the  design  of 
dynamic  compensators  is  in  general  reduced  to  the  selection 
of  pairwise-orthogonal  eigenvectors  and  reciprocal  eigenvectors 
for  two  families  of  well-defined  subspaces  which  are  para- 
metrised by  associated  self-conjugate  eigenvalue  spectra. 

This  selection  can  be  effected  by  the  use  of  a powerful  new 
algorithm^^  (see  Appendices  10  and  11)  which  requires 

- — — - — 
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the  performance  of  restricted  elementary  row  and  column 
operations  on  matrices  formed  from  the  spanning  vectors  of 
these  subspaces.  The  principal  computational  attraction  of 
this  algorithm  is  that  no  operations  with  polynomial  matrices 
are  involved,  so  that  dynamic  compensators  for  large-scale 
systems  can  be  readily  designed. 


3.  SOFTWARE  PACKAGE  DEVELOPMENT 


3.1  Capability  of  EIGENFORTRAC  Software 


The  EIGENFORTRAC  software  package  is  essentially  an  up- 

(18) 

dated  version  of  FORTRAC  ' based  solely  on  the  unifying 
method  of  entire  eigenstructure  assignment  and,  in  particular, 
on  the  powerful  algorithm  for  the  design  of  dynamic  com- 

(16)  (17)  P • ■ MB  4»  In  ^ M 4 M 4b  A a Vn  M 4 * « A A ^ A A 4b  4b  A A A JVn  an  a!# 


pensators 


. Synthesis  techniques  for  state-feedback 


regulators,  observers,  output- feedback  regulators,  and 
dynamic  compensators  are  embodied  in  EIGENFORTRAC.  These 
techniques  have  been  applied  to  the  design  of  controllers 
for  a variety  of  aircraft  in  a number  of  flight  modes.  Thus, 
for  example,  digital  flight  control  systems  have  been 
designed  by  D'Azzo  and  Porter for  the  F-4  fighter  air- 
craft and  by  D'Azzo  and  Kennedy^20^  for  the  C-141  transport 
aircraft. 


3.2  Configuration  of  EIGENFORTRAC  Software 


The  EIGENFORTRAC  program  configuration  has  been  described 
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by  Porter,  Bradshaw,  and  Daintith * ^ , together  with  a 
detailed  description  of  all  the  EIGENFORTRAC  subroutines. 
Detailed  listings  of  the  computer  output  for  a simple 
example  illustrating  the  design  of  discrete-time  tracking 
systems  incorporating  error-actuated  dynamic  compensators 
have  also  been  provided 

3. 3 Operation  of  EIGENFORTRAC  Software 

The  basic  requirements  of  EIGENFORTRAC  are  the  plant,  input, 
and  output  matrices  (A,B,C)  of  the  uncontrolled  system 
described  in  the  continuous-time  domain.  The  class  of 
controller  required  is  then  specified  (eg,  state-feedback 
regulator,  output- feedback  regulator,  error-actuated  dynamic 
compensator)  together  with  the  sampling  interval,  T.  The 
plant,  input,  and  output  matrices  (A(T) ,B(T) ,C)  of  the 
sampled  uncontrolled  system  are  then  computed,  and  appro- 
priate augmentation^  is  automatically  introduced.  The 
closed-loop  eigenvector  and  reciprocal  eigenvector  subspaces 
are  then  computed,  and  pairwise-orthogonal  sets  of  closed- 
loop  eigenvectors  and  reciprocal  eigenvectors  are  then 
selected  from  these  subspaces.  Finally,  these  sets  of 
eigenvectors  and  reciprocal  eigenvectors  are  used  in  the 
computation  of  the  compensator  matrices.  The  performance  of 
the  resulting  controller  is  checked  by  performing  a discrete- 
time simulation  which  is  followed  by  a continuous -time 
simulation  using  a Kutta-Merson  routine  in  the  case  of 
promising  designs. 
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4.  CONCLUSIONS 

Fundamental  new  insights  into  the  structure  of  linear  multi- 
variable  systems  have  been  obtained  by  developing  a unified 
theory  of  entire  eigenstructure  assignment.  These  system- 
theoretic  results  have  been  implemented  in  the  comprehensive 
software  package  E I GEN FORT RAC ^ which  is  currently  available 
for  the  routine  computer-aided  design  of  on-board  digital 
controllers  for  a wide  range  of  flight-control  applications. 
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Abstract 

It  is  shown  that  the  synthesis  of  closed-loop 
linear  multivariable  discrete-time  systems  can  be 
directly  effected  by  performing  equivalence  transfor- 
mations on  appropriate  polynomial  matrices.  These  poly- 
nomial matrices  are  the  Smith  canonical  forms  of  the 
closed-loop  characteristic  matrices  of  such  systems  sub- 
ject to  the  constraints  imposed  by  the  fundamental 
theorem  of  linear  state-variable  feedback. 
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1,  Introduction 


In  this  paper  it  is  shown  that  the  synthesis  of 
closed-loop  linear  multivariable  discrete-time  systems 
governed  by  state  and  feedback  equations  of  the  respective 
forms 

x(k+l)  = A::(k)  + Bu(k)  (la) 

* •# 

and 

u(k)  - Gx(k)  (lb) 

where  x(k)£  Rn,  u(k)£  Rm,  and  (A,B)  is  a reachable  pair 
can  be  directly  effected  by  performing  equivalence  trans- 
formations on  appropriate  polynomial  matrices.  This 
synthesis  procedure  consists  essentially  in  the  specifica- 
tion of  the  Smith  canonical  forn^  S(A)  of  the  closed- 
loop  characteristic  matrix  (AI^-A-BG)  subject  to  the 
constraints  on  S(A)  imposed  by  the  fundamental  theorem  of 
linear  state-variable  feedback ^ ^ . 

2,  Synthesis  Procedure 

Thu3,  if 

S ( A ) = diag  (1»1»....»  ( A ) » (A)  * • • • . f ^(M  # ( A ) ) 

(2) 

where  the  ^(A)  (i=l,2, . . . . ,q)  are  any  monic  polynomials 
in  R[a]  such  that 
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(3) 


and 

q 

E deg^ . (X)  - n , (4) 

i-1  1 

then  there  exists  a matrix  g£  so  that  S(X)  is  the 

Smith  canonical  form  of  (XI  -A-BG)  orovided  that 

~ n - - - 


P p 

Z degij/ . (A)  > Z <,  (l<P<q)  (5) 

i-1  1 “i-1  1 

where  *^,><2^  are  the  ordered  Kronecker  invariants'”^ 


of  the  pair  (A,B) . It  is  accordingly  evident  that  closed- 
loop  discreto-time  systems  governed  by  state  and  feedback 
equations  of  the  form  (1)  can  be  synthesised  by  the 
following  procedure: 

(i)  Determine  the  Kronecker  invariants  <i (i-1, 2, . . . . ,m) 

of  the  pair  (A#B); 

-»  **» 

(ii)  Prescribe  an  admissible  matrix  S(X)  on  the 
basis  of  the  values  of  the  <i (i-1, 2, . . . . ,m) ; 

(iii)  Transform  S(X)  by  elementary  row  and  column 

operations  into  an  equivalent  polynomial  matrix 
of  the  form 

Z ( X)  - XI  -A-BG  j (6) 

(iv)  Determine  the  set  of  linear  simultaneous 
equations  satisfied  by  the  elements  of  the 
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feedback  matrix  G by  inspection  of  E(X); 

(v)  Solve  the  set  of  linear  simultaneous  equations 
for  the  elements  of  the  feedback  matrix  G. 


3.  Illustrative  Example 


This  procedure  can  be  conveniently  illustrated  by 
synthesising  a closed-loop  system  governed  by  the  state 
and  feedback  equations 


x(k+l)  “ 

'O  , 1 , 2 

-2  , 3,0 

x (k)  + 

‘l  , 2 

1 , o 

-2  ,-1,0. 

O , 0. 

u ( k ) 


(7a) 


and 


u (k) 


Gx  (k) 


gll 

/ s12 

' g13 

•g21 

t 922 

' g23 

x(k) 


(7b) 


so  that  the  eigenvalues  of  the  closed-loop  plant  matrix 
are  all  equal  to  zero.  In  this  case  it  is  evident  from 
equations  (7)  that 


E(X) 


X~gll“2g21  ' “1_g12'2g22  ' “2”g13"2g23 
2“gH  r *“3-g12  , -gL3 

2 , 1 , X 


(3) 


and  that  the  associated  Kronecker  invariants  are  <^*2, 
<2»1.  The  conditions  (3),  (4),  and  (5)  therefore  indicate 


that 
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(X)  • diag(l,X,\2)  (9) 

and 

S2(X)  - diag(l,l,X3)  (10) 

are  the  only  admissible  forms  of  the  Smith  canonical  form 
S(X)  of  the  characteristic  matrix  of  the  closed-loop  system 
governed  by  equations  (7). 

In  the  case  of  S^(X)  it  is  readily  found  that 


0 , 1 , o' 

1,0,  o' 

•2  , 1 , X] 

X , -2  , -1 

0 , X , 0 

o 

* 

o 

♦ 

1 

H* 

* 

o 

0 

1  

1 

<M 

*< 

O 

o 

1 

I 

r- 1 

% 

o 

* 

_2i 

‘x  , o , o' 

0 , X , 0 

- ^(XMll) 

.2  , 1 , X< 

so  that  comparison  of  equations  (8)  and  (11)  indicates 
that  the  corresponding  feedback  matrix  in  equation  (7b) 
is 


\ 


2 , -3  , 0 


(12) 


Similarly,  in  the  case  of  S2(X),  it  is  readily  found 
that 


U 

0 

0 
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'1  , 0 , cl  p.  * o , ol  ~ x ,0,1 

“2  , X , 1 0,1,0  2 - X ~ ,1,0 


x ,1,01  Lo  , 0 , X 


1 r 0 , 01 


x , 0 , i' 
0 , X , —2 
.2,1,  X. 


t2(\)  (13) 


so  that  comparison  of  equations  (8)  and  (13)  indicates 
that  the  corresponding  feedback  matrix  in  equation  (7b) 


2,-3,  2 


-1»  1 / -5/2 


It  is  clear  that,  as  desired,  the  characteristic  polynomial 
of  the  closed-loop  plant  matrix  is 


c(X) 


in  both  cases  but  that  the  minimum  polynomials  associated 
with  the  feedback  matrices  G^  and  G2  are  respectively 

2 

mx(X)  = X (16) 


m2  (X) 


4.  Conclusion 


This  procedure  for  the  synthesis  of  linear  multivariable 
discrete-time  feedback  systems  constitutes  a generalised 


I 


- 17  - 

L/ 

eigenvalue-assignment  procedure  in  t ha t.  both  the  cyclic 
structure  and  the  eigenvalues  of  the  closed-icop  plant 
matrices  are  synthesised.  Moreover,  the  fact  that  the 
synthesis  of  such  systems  is  directly  effected  by  performing 
equivalence  transformations  on  S(,\)  ensures  that  only  those 
cyclic  structures  which  are  conformable  with  the  constraints 
imposed  by  the  fundamental  theorem  of  linear  state-variable 
feedback  are  considered.  In  particular,  tha  synthesis 
procedure  facilitates  the  assignment  of  both  closed-loop 
characteristic  polynomials  and  admissible  closed-loop 
minimum  polynomials.  This  facility  is  particularly  important 
in  the  case  of  discrete-time  systems  since  it  obviously 
provides  a basis  for  the  design  of  time-optimal  linear 

(4\ 

multivariable  control  systems  . It  is  evident,  however, 
that  the  generalised  eigenvalue-assignment  procedure  is 
equally  applicable  to  the  synthesis  of  linear  multi1 variable 
continuous-time  feedback  systems. 
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Abstract 

The  algorithm  for  the  computation  of  a basis  for 
ker  [A“*0In  , B]  presented  in  this  paper  greatly  facilitates 
the  synthesis  of  state-feedback  regulators  by  entire  eigen- 
structure  assignment.  It  is  ideally  suited  for  digital 
computer  implementation  and  can  be  readily  dualised  for  use 
in  the  synthesis  of  full-order  observers  by  entire  eigen- 
structure  assignment. 


I 


l 


21 


u 


1.  Introduction 


In  view  of  the  recent  results  obtained  by  Kimura^ 

(2) 

and  Moore  , it  is  evident  that  an  efficient  algorithm  for 
the  computation  of  a basis  for 

ker  S ( X ) = ker  [a-X  I , b‘|  (!) 

-*  o — o— n — 

where  Xq£  C and  [A-XQIn  , b}  £ ^nx(n+m)  essential  for 
the  synthesis  by  entire  eigenstructure  assignment  of  state- 
feedback  regulators  for  multivariable  linear  systems  governed 
by  state,  output,  and  control-law  equations  of  the  respective 
forms 


* 1 


; 

n 

P 

i 

| 


x(t) 

= Ax  ( t)  + Bu(t)  , 

(2) 

y(t) 

= Cx  ( t)  , 

(3) 

and 

u(t) 

= Kx(t) 

(4) 

Indeed,  the 

real  state- feedback 

matrix(1)(2> 

K = 

[«1  • ' * * * ' 

[Xi  , X2  f ...  , *n 

l"1  (5) 

simultaneously  assigns  the  self-conjugate  distinct  eigen- 
value spectrum  {X^,Xj, . . . , Xr}  and  the  corresponding  eigen- 
vector set  {X^,  X'2 1 • • • / Xn)  to  the  closed-loop  plant  matrix 
(A+BK) £ Rnxn  just  in  case 


1 

; 


[*i  ' “i]'G  ker  (i*l, 2, . . . ,n) 
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2.  Algorithm 


The  following  algorithm  provides  an  efficient  means 

for  the  computation  of  a basis  for  ker  [A-X  I , d] s 

*•»  o ~ n •» 

(i)  Form  the  matrix 


!<V  - 


A-X  I„  , B 
- o-n 

Jn+m 


where  X^  £ C ; 


(ii)  Perform  elementary  column  operations  on  S(X  ) until 

- o 


su j . : 

- o ~ 


?u  ' 9 


?21  ' f2 


S(XJ 

"»  O 


where  Cnxr,  rank  « r ■»  rank  [a— AoT.n  , b]  , 

i21E  C(n+m)xr,  and  S22E  C{n+tn)x(n+m"r^. 

The  required  basis  vectors  for  ker  [a-X  I , d]  are  then 

~ o~n  — 

given  by  the  (n+m-r)  columns  of  S22#  where  obviously  r ® n 
in  case  XQ  is  not  an  input-decoupling  zero  of  the  system. 
This  follows  from  the  fact  that  equations  (7)  and  (8)  imply 
that 


S(X0)  - S(Xo)E  • S(Xo)  [5l  , E2] 


^"XoJn  ' ?3fl  • 


[A-X  1 , 13]E. 

- o-n  - J •» . 


where  E^  ^(n+m)x(n+m)  a pr0(juct  Qf  elementary  matrices, 
rank  - r,  rank  E2  » n + m - r,  [A-XQin  , - Sn, 

^‘Xo*n  ' ?-^?2  " 9*  5l  " ?21'  and  ?2  “ ?22* 

3.  Illustrative  Example 


This  algorithm  can  be  conveniently  illustrated  by 
synthesising  a state-feedback  regulator  for  a multivariable 
linear  system  characterised  by  the  matrices 


-1»0,0' 
0,-2,  0 
.1  , 1 , -3 


(10) 


B 


1 , 0 
0 , 1 
IP  , 1 


(11) 


and 


C - 


1,2,0 

1,2,1 


(12) 


which  is  such  that 


O (A+BK)  - (-1,-2, -3}  - 0(A) 


(13) 


but  such  that  the  'slow'  mode  corresponding  to  the  eigen- 
value ■ -1  is  eliminated  from  the  output.  Hence,  in 
accordance  with  the  algorithm,  it  is  found  that 


25 


I 

In  view  of  the  equivalences  (14),  (15),  and  (16)  it  therefore 
follows  from  the  algorithm  that 


ker  S (-1) 


ker  S(-2) 


span  ( 


2 

0 

1 - 
0 
0 


» 


span  s 


» 


(17) 


(18) 


and 


ker  S(-3) 


(19) 


It  is  evident  from  ker  S(-l),  ker  S(-2)  , and  ker  S(-3)  that 
the  closed- loop  eigenvectors  corresponding  to  the  eigenvalue 
spectrum  {-1,-2, -3}  * c(A+BK)  can  be  assigned  to  the  respective 
subspaces 


E(-l) 


span 


i 


(20) 


subject  only  to  the  requirement  that  the  resulting  set  of 


eigenvectors  be  linearly  independent.  In  case 

l*{  / «J]'  * [-2  , 1 , 0 , 0 , 1]'  , (23) 

l>2  ' “Jl'  = [l  , 0 , 1 , -1  , 0]'  , (24) 

and 

[XJ  , wj]"  - [0  , 1 , 0 , O , -1]'  , (25) 

it  follows  from  equation  (5)  that 


(26) 

and  therefore  from  equations  (10) , (11),  and  (2G)  that 


The  eigenvalues  have  accordingly  been  unaltered  by  state 
feedback/  as  required,  but  the  corresponding  eigenvectors 
have  become 


{x1fX2»x3} 


r 

’-2 

T 

0 

1 

9 

0 

9 

1 

> 

_0_ 

.1 

p 

where 


Xx  = [-2,1,  0]'£ker  C’  , 


(28) 


(29) 


as  required.  This  elimination  of  the  'slow'  mode  correspond- 
ing to  the  eigenvalue  = -1  is  possible  because  -1  is  an 

invariant  zero  and  X-^  is  a corresponding  state  zero-direction 
(3) 

of  the  system 


4.  Conclusion 

This  algorithm  for  the  computation  of  a basis  for 
ker  [A-A  In  , B]  greatly  facilitates  the  synthesis  of  state- 
feedback  regulators  by  entire  eigenstructure  assignment  since 
it  is  ideally  suited  to  digital  computer  implementation.  In 
addition,  it  is  evident  that  the  same  algorithm  also  greatly 
facilitates  the  synthesis  of  full-order  observers  by  entire 
eigenstructure  assignment  since  it  can  clearly  be  used  for 
the  computation  of  a basis  for  ker  [A'-A^^  , C'] . 
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Abstract 

In  this  paper,  results  are  presented  which  facilitate 
the  complete  exploitation  cf  state  feedback  in  the  assignment 
of  the  entire  closed-loop  eigenstructure  of  multivariable 
linear  systems.  These  results  include  an  algorithm  for  the 
direct  computation  of  the  state- feedback  matrix  which  assigns 
prescribed  Jordan  canonical  forms,  eigenvectors,  and  generalised 
eigenvectors  to  the  plant  matrices  of  closed-loop  systems. 

This  algorithm  is  illustrated  by  assigning  the  entire  closed- 
loop  eigenstructure  of  a third-order  two-input  discrete-time 
system  in  such  a way  that  the  resulting  closed-loop  system 
exhibits  finite  settling  time  behaviour. 


It  is  well  known  that,  except  in  the  case  of  single- 
input systems,  specification  of  closed-loop  eigenvalues  does 
not  define  a unique  closed-loop  system.  This  non-uniqueness 
has,  however,  been  only  partially  exploited  in  only  a few 
instances  by  algorithms  which  permit  the  specification  of  a 
number  of  components  of  the  closed-loop  eigenvectors  (Srinathkumar 
and  Rhoten  1975,  Shah  et  al  1975)  and  by  algorithms  which  avoid 
large  feedback  gains  (Porter  and  Crossley  1972,  Lee  1975).  The 
results  presented  in  this  paper  facilitate  the  complete  exploit- 
ation of  state  feedback  in  the  assignment  of  the  entire  closed- 
loop  eigenstructure  of  multivariable  linear  systems.  These 
results  include  an  algorithm  for  the  direct  computation  of  the 
state-feedback  matrix  which  assigns  prescribed  Jordan  canonical 
forms,  eigenvectors,  and  generalised  eigenvectors  to  the  plant 
matrices  of  closed-loop  systems.  The  expression  for  this  state- 
feedback  matrix  assumes  a simple  form  which  is  equivalent  to 
that  obtained  by  Kimura  (1975)  in  the  context  of  eigenvalue 
assignment  by  output  feedback  and  by  Mo6re  (1976)  in  the  context 
of  output  regulation  for  the  special  case  of  distinct  prescribed 
eigenvalues . 


t 


i 


l 


1 


[a-x1i  , b] 


Yx 


Yx 
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(2,  j) 
i 

(2,j) 


L t 


Yx 


(1,3) 


[A-X1I  , B] 


(m.  . , j)' 


YX' 


w <mji'j) 

L'Ai 


'•Aj 


(lb) 


(1V 


I l. 


L 


I 


( j*l , 2 , • • . / k^  ; i— l/2,...,p) 

generate  k^  strings  of  vectors  associated  with  the  eigenvalue 

X , where  v.  ,J/  is  the  £th  vector  in  the  jth  string  of 
i 

length  associated  with  the  eigenvalue  X^.  The  vectors 
Yx*1'^  ( j*l,2, . . . are  the  eigenvectors  associated 
with  the  eigenvalue  X^,  whilst  the  remaining  vectors  in  each 
of  the  k^  strings  of  vectors  generated  by  equations  (1)  are 
generalised  eigenvectors  associated  with  the  eigenvalue  X^. 
The  total  number  of  vectors  associated  with  the  eigenvalue  X^ 
is  evidently 


m. 


i 
E m 
j-1 


ji 


(i=l,2, . . . ,p) 


(2) 


and  the  entire  set  of  vectors  associated  with  the  eigenvalue 
spectrum  {X^, » • • • * X^}  will  accordingly  serve  as  a basis 
for  n-dimensional  state  space  only  if 


P 

E m. 
i-1  1 


(3) 


In  case  the  eigenvalues  A^  and  the  integers  and  k^  are 

chosen  so  that  this  entire  set  of  vectors  not  only  satisfies 
(3)  but  is  also  linearly  independent  and  self-conjugate,  then 
the  real  state-feedback  matrix 

k=  [w/1'11 w/Vp'V]  iV1'15 vx(V'kp)] 

- ' A1  ~ Ap  P ' Ai  ' Ap  P 


is  such  that  the  Jordan  canonical  form  of  the  nxn  closed-loop 
plant  matrix  (A+BK)  contains  the  eigenvalue  A i ( i=l , 2 , . . . ,p) 
with  geometric  multiplicity  k^  and  algebraic  multiplicity  m^. 

This  follows  from  the  fact  that  if  the  real  state-feedback 
matrix  K is  such  that  the  Jordan  canonical  form  of  the  closed- 
loop  plant  matrix  (A+BK)  contains  the  eigenvalue  A^^  (i=l , 2 , . . . ,p) 
with  geometric  multiplicity  k^  algebraic  multiplicity  nu  , 
and  associated  eigenvectors  ( j=l, 2 , . . . ,k^)  then 


[A- A^  , B] 


(l,j) 


(1,  j) 


“ - 0 


&-xi*  ' 5.1  ' i 


(2 , j ) 


(2 , j ) 


(1  r 3) 


[A-Xjl  , B]  1 


*•  A j 


Kv.  '"'ji'31 

A ^ * 


v . (n'j  i~ 1 f j ) , 

~Ai 


(51^) 


0 


34 


■ 


p* 


where  the  m^  satisfy  (2)  and  (3) . 

It  is  evident  that,  in  the  special  case  when  p * n and 
ki  ■ n»i  ■ 1 (i*l,2, . . . ,n) , then  j « 1 and  m ^ i * 1 (i*l, 2, . . . ,n) . 
Each  of  the  sequences  of  equations  (1)  accordingly  reduces 
to  just  a single  equation,  and  there  are  clearly  n such 
equations 


v (1,D 

~ i 

w (1'1J 


0 


( i=l, 2 , . . . , n)  (6) 


for  the  eigenvectors  v,  ' (i=l , 2 , . . . ,n)  of  the  nxn  closed- 

i 

loop  plant  matrix  (A+BK) . In  thxs  special  case,  the 
expression  (4)  for  the  state-feedback  matrix  assumes  the 
simple  form 


k - ktl'1,.«x(l'11 (i.inr  (un,  a,D (i.i)i 

-*2  "An  J L'*l  -12  -xn  j 


-1 


(7) 


which  is  equivalent  to  that  obtained  by  Kimura  (1975)  and 
Moore  (1976) . The  computation  of  K in  the  case  of  distinct 
eigenvalues  thus  reduces  to  the  determination  of  the  kernels 
of  each  of  the  n matrices. 


?A 


i 


(i=l, 2 , . . . ,n)  . (8) 


3,  Illustrative  Example 

These  results  can  be  conveniently  illustrated  by  assigning 
the  entire  closed- loop  eigenstructure  of  the  discrete-time 
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T 


system  governed  by  the  state  and  feedback  equations  (Porter  1976) 


f* 

x(k+l)  - 

"0  , 1 ,2' 

-2,3,0 

x(k)  + 

1 , 2‘ 

1 , o 

r 

1 

0 

•H 

1 

CN 

1 i 

1 

O 

Qi 

u(k) 


and 

i 

u(k)  - Kx (k) 


(9a) 


(9b) 


V 


l 


in  such  a way  that  the  resulting  closed-loop  system  exhibits 
finite  settling  time  behaviour.  Such  an  assignment  clearly 
requires  that  p « 1,  m^  ■ 3,  and  “ 0 but,  in  consonance 
with  the  fundamental  theorem  of  linear  state-variable  feedback 
(Dickinson  1974) , it  is  possible  further  to  require  either 
that  k^  » 2,  m^  * 2,  and  m21  - 1 or  that  k^  * 1 and  m^  = 3. 
In  the  former  case,  equations  (1)  indicate  that 


(10) 


I 


constitutes  an  admissible  set  of  closed-loop  eigenvectors  and 
generalised  eigenvectors  whilst,  in  the  latter  case,  equations 
(1)  indicate  that 


<v  <l'« 
"A1 


12,1) 

~A1 


I 


(11) 


36 


constitutes  an  admissible  set  of  closed-loop  eigenvectors  and 
generalised  eigenvectors  since  also 


'0  , 1 , 2 , 1 , 2 
-2  , 3 ,0,1,0 

-2  , -1  ,0,0,0. 


and  therefore 


Ker 


* span 


In  the  former  case, 


(1,1)  (2,1.)  w <1,2>} 

, w.  , w,  j 


<WX  * WX 

"A1 


ta  ■ p ■ m 


so  that  in  view  of  equation  (4) 


0 , 11  , 8 

-1  * "6  , -3. 


0,1,1' 
0,-3,  -2 
1,2,0 


and  therefore 


0 , 0 , 0 1 


( A+BK)  - 0,0,0 

-2  ,-1,0 


2 , -3  , 0 

-1,1,  -1 
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i li 


which  has  the  Jordan  canonical  form 


0,1,0 


0,0,0 

0,0,0 


together  with  the  eigenvectors  and  v.^'2^  and  the 

(2  1)  1 1 

generalised  eigenvector  v^ ' ' ' prescribed  in  equation  (10), 


a3  required.  In  the  latter  case, 


(1,1)  (2,1)  (3,1), 

, w\  » w\  ' 

1 ~A1  ~A1 


■ m 


so  that  in  view  of  equation  (4) 


which  has  the  Jordan  canonical  form 

fo  , 1 , 01 


- 38  - 

eigenvectors  v^^'^  and  v^^'^  prescribed  in  equation  (11), 
as  required. 

4.  Conclusion 

These  results  facilitate  the  complete  exploitation  of 
state  feedback  in  the  assignment  of  the  entire  closed-loop 
eigenstructure  of  multivariable  linear  systems  and  are  clearly 
equally  applicable  to  both  continuous-time  and  discrete-time 
systems.  It  is  evident  that,  even  in  the  case  of  systems  for 
which  the  pair  (A,B)  is  uncontrollable,  certain  prescribed 
eigenvectors  of  (A+BK)  can  be  assigned  by  state  feedback.  In 
the  case  of  systems  with  asymptotically  stable  but  uncontrollable 
modes,  it  is  therefore  frequently  possible  to  achieve  significant 
improvements  in  the  dynamical  behaviour  of  such  systems  by  the 
introduction  of  appropriate  state-feedback  controllers. 
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Abstract 

In  this  paper,  an  algorithm  is  presented  which  greatly 
facilitates  the  complete  exploitation  of  state  feedback  in 
the  assignment  of  the  entire  closed-loop  eigenstructure  of 
controllable  multi-input  systems.  This  algorithm  is  a 
generalisation  of  the  algorithm  of  MacLane  and  Birkhoff  (1968) 
for  the  computation  of  a basis  for  the  null  space  of  a matrix 
and  is  ideally  suited  to  digital  computer  implementation. 

The  algorithm  readily  yields  the  vectors  which  are  required 
(Porter  and  D'Azzo  1977)  for  the  simultaneous  assignment  of 
Jordan  canonical  forms,  eigenvectors,  and  generalised  eigen- 
vectors to  the  plant  matrices  of  closed-loop  controllable 
multivariable  linear  systems.  The  effectiveness  of  the 
algorithm  is  illustrated  by  assigning  the  entire  closed-loop 
eigenstructure  of  a third-order  two-input  discrete-time  system 
in  such  a way  that  the  resulting  closed- loop  system  exhibits 
time-optimal  behaviour. 
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1.  Introduction 


The  algorithm  presented  in  this  paper  readily  yields 
the  vectors  which  are  required  (Porter  and  D'Azzo  1977) 
for  the  determination  of  the  state- feedback  matrix  which 
simultaneously  assigns  Jordan  canonical  forms,  eigenvectors, 
and  generalised  eigenvectors  to  the  plant  matrices  of  closed- 
loop  controllable  multi-input  linear  systems.  These  vectors 
satisfy  the  sequences  of  equations  (Porter  and  D'Azzo  1977) 


[A-XiIn  , B] 


rv. (1'j) 


(i»  j) 
u’"i 


Yx 


rvx(2'j) 

~XL 


w 


Yx 


( l * j ) 


(la) 


(lb) 


In  , §] 


Yx 


Yx 


(mjt, j) 


(mjl-l,j) 


(lm^) 


( j “1 , 2 , . . . , f i“l , 2 , . • . , p) 

which  together  generate  k^  strings  of  eigenvectors  and 

generalised  eigenvectors  associated  with  the*  eigenvalue  X^ 

(l  i ) 

where  v.  ,J  is  the  tth  vector  in  the  jth  string  of  length 
i 

m^  associated  with  the  eigenvalue  X^.  In  case  the  eigenvalues 
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X^  (i-1,2, . . . ,p)  and  the  integers  m^  and  k^  are  chosen  so 
that  this  entire  set  of  eigenvectors  and  generalised  eigen- 
vectors is  linearly  independent  and  self-conjugate,  then  the 
real  state-feedback  matrix  (Porter  and  D'Azzo  1977) 


(m.  ,k) 

V P 1 L (l.D 


r ‘V'Vi-i 

P 


is  such  that  the  Jordan  canonical  form  of  the  nxn  closed- 

loop  plant  matrix  (A+BK)  contains  the  eigenvalue  X.  (i«l , 2 , . . . ,p) 

with  geometric  multiplicity  and  algebraic  multiplicity 


( i“l , 2 , . . . , p) 


It  is  evident  that,  in  the  special  case  when  p = n and 


* m^  * 1 (i=l, 2 , . . . ,n) , then  j * 1 , 


1 (i*l, 2, . . . ,n) , 


and  each  of  the  sequences  of  equations  (1)  reduces  to  just 

a single  equation  for  the  eigenvector  v of  the  nxn 

i 

closed-loop  plant  matrix  (A+BK)  associated  with  the  eigen- 
value  X^.  In  this  special  case  of  self-conjugate  distinct 
eigenvalues  {X^, X^, . . . , X^} , the  general  expression  (2)  for 
the  state-feedback  matrix  K assumes  the  simple  form  obtained 
by  Kimura  (1975)  and  Moore  (1976) . 


2,  Algorithm 


<"*  .p'V 

v.  p } and 
~AP 


The  vectors  {v, 
~*1 


9 • • • 9 
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(mk  P'V 

(1.1) 

fy.  , ...«w  } required  for  the  determination  of 

~A1  'Ap 

the  state-feedback  matrix  K expressed  by  equation  (2)  can 
be  readily  computed  by  the  following  algorithm  which  is  a 
generalisation  of  the  algorithm  of  MacLane  and  Birkhoff 
(1968)  for  the  computation  of  a basis  for  the  null-space 
of  a matrix: 

(i)  Form  the  matrix 


?<v 


?n+m 


b1 


(4) 


for  \ ( i= 1 » 2 , . , . , p)  ; 


(ii)  Perform  elementary  column  operations  on  S(\^)  until 


S(*i> 


fc  (1.3)  o (l.j)- 
'-11  ' §12 


q (1.3)  q (1.3) 
§21  ' §22 


(5) 


where  S^1'^  Cnw\ 
n * rank  |^A— , b] 
pair; 


§121,j)  = °'  and  rank  - 
since  (A, B)  is  a controllable 

«N» 


(iii)  Perform  successive  elementary  column  operations  on 

S(1'j) (A.)  until 
- 1 


* (m. .-1, j) 

S J1 


fu1'3’  • Yx'1"31'1'3’1 


S (l.j)  s (mji'^ 
.-21  ' §22 


s(mji'j) 


(6) 


where 


Yx 


wXj 


i 

(m^O) 


(mji.J) 

-22 


(7) 


(»4  4 * j) 


The  matrices  (V,  VJ‘’  Jl  , V,  v*  ' Jl  , . . . , VA  ^ } and 


d,j>  v (2,j) 

1 i 

{wx(1,j)*wx(2,j) W,  ^ } thus  generated  are  clearlv 

i ~Ai  ~Ai 

such  that 


[A-^In  , B] 


rv  UO> 
~Ai 


w, 


(lo) 


(8a) 


(A-Xiln  , BJ 


rv  <2*:n 

~Xi 


wx(2'j) 

■ i 


Yx 


do) 


(8b) 


[A-X1In  r B] 


V,  J1  1 
~Ai 


L?X 


(m.  -1, j) 

Yx,  31 


(8m^) 


and  are  therefore  such  that  the  vectors 


* ( 

L-xi 


(njt,j)'  (n^O) 


» w, 


j i 1 


(Hji*1, 2'  ‘ * * required  for  use  in  equation  (2)  are  linear 


ji 


combinations  of  corresponding  columns  of  successive  members 

r (n . , , j)  ' (n., f j) S , 

of  the  entire  sequence  of  matrices  |V^  J 7 


(n^-1,2, 


' ^ j i^  * 
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3.  Illustrative  Example 

This  algorithm  can  be  conveniently  illustrated  by 
assigning  the  entire  closed-loop  eigenstructure  of  the  discrete- 
time system  governed  by  the  state  and  feedback  equations 
(Porter  1976a, b) 
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and 

'l  / 2 , 1 , 0 r O' 

3/0,1#  0,0 

"1  ,0,0,  0 , 0 

0,0,0,  0,  1 
1 ,0,0,  0 , -2 
0,1,0,  1 , 0 
0,0,1,  0 , 8 
0 , 0 , 0 , -1  , -3 


1 , 2 , 1 , 0 , 1 
3 ,0,1,  0 , -2 
“1  ,0,0,  1 , 0 


0,0,0,  1 , 0 
1 , 0 , 0 , -3  , 0 
0,1,0,  2 ,3/2 

0 , 0 , 1 , 11  , -2 

0 ,0,0,  *6  , 0 


(10b) 


It  is  evident  from  the  equivalences  (10a)  and  (10b)  that 

~A1  ~A1  ~A1 

constitutes  an  admissible  set  of  closed-loop  eigenvectors 
and  generalised  eigenvectors  and,  correspondingly,  that 

« 

■.  (12) 

. 


"X, 


0 

n 

8 

i 

/ 

-1 

-6 

-3 

-3 


1 

-2 

0 


(ID 


The  required  state-feedback  matrix  determined  by  equation 
(2)  is  therefore 


U , 2 , 0 


(13) 
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so  that  the  plant  matrix  of  the  closed-loop  system  governed 
by  equations  (9)  and  (13)  is 


A + BK 


0 , 0,0“ 
0,0,0 


L-2  , -1  , Oj 


(14) 


This  plant  matrix  is  clearly  nilpotent  of  index  two,  as. 
required,  and  furthermore  possesses  the  eigenvectors  and 
generalised  eigenvectors  prescribed  in  equation  (11) . 


4.  Conclusion 


In  this  paper,  an  algorithm  has  been  presented  which 
greatly  facilitates  the  synthesis  of  state-feedback  regulators 
by  entire  eigenstructure  assignment.  This  algorithm,  which 
is  equally  applicable  to  both  continuous-time  and  discrete- 
time systems,  has  been  illustrated  by  assigning  the  entire 
closed-loop  eigenstructure  of  a third-order  two-input  discrete- 
time system  in  such  a way  that  the  closed-loop  system  exhibits 
time-optimal  behaviour.  In  view  of  the  simple  elementary 
column  operations  involved,  it  is  evident  that  the  algorithm 
is  ideally  suited  to  digital  computer  implementation. 
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Abstract 

In  this  paper,  the  method  of  entire  eigenstructure 
assignment  (Kimura  1975,  Moore  1976,  Porter  and  D'Azzo  1977) 
is  applied  to  the  design  of  linear  multivariable  continuous- 
time output-feedback  regulators.  It  is  shown  that,  in  the 
case  of  self-conjugate  distinct  eigenvalue  spectra,  the  closed- 
loop  eigenstructure  assignable  by  output  feedback  is  constrained 
by  the  requirement  that  the  eigenvectors  and  reciprocal  eigen- 
vectors lie  in  well-defined  subspaces.  The  method  is  illustrated 
by  designing  an  output- feedback  regulator  for  a third-order 
continuous-time  system. 
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1.  Introduction 


In  most  practical  cases,  it  is  impossible  to  implement 
state-feedback  control  laws  since  the  state  of  the  plant  is 
inaccessible  and  only  the  plant  output  is  available  for 
control  purposes.  Much  effort  (see,  for  example,  Davison 
and  Wang  (1975))  has  accordingly  been  expended  on  the  investi- 
gation of  the  closed-loop  dynamics  achievable  by  the  implement- 
ation of  output- feedback  control  laws.  However,  apart  from 
the  partial  results  obtained  by  Kimura  (1975) , this  effort 
has  led  to  results  concerned  only  with  closed-loop  eigen- 
values and  not  with  closed-loop  eigenvectors.  In  this  paper, 
the  method  of  entire  eigenstructure  assignment  (Kimura  1975, 
Moore  1976,  Porter  and  D'Azzo  1977)  is  therefore  applied  to 
the  design  of  output- feedback  regulators  for  multivariable 
linear  continuous-time  systems  governed  by  state  and  output 
equations  of  the  respective  forms 

X ( t ) - Ax ( t)  + Bu(t)  (1) 

+0  «•  <*»  «w 

and 

^(t)  - Cx(t)  , (2) 

where  A£lRnxn,  b£  Rnxm,  R*5*11,  rank  B = m,  and  rank 

C * p. 

2,  Theory 

Thus,  if  output  feedback  is  applied  to  the  system 
governed  by  the  state  equation  (1)  in  accordance  with  the 


\ 

I 


I 

[ 


control-law  equation 


u(t)  - Gy(t)  (3) 

and  the  output- feedback  matrix  G £ Rmxp  is  such  that  the 
closed-loop  plant  matrix  (A-»BGC)  has  a self-con  jugate 
spectrum  {X^, Xj# • • • # Xn)  of  distinct  eigenvalues  and  correspond 
ing  eigenvector  and  reciprocal  eigenvector  sets  (Xi *Xo * • • • » x } 
and  {^,<(>2# . . . , $n},  then  obviously 

(-“Xii+??S) *i  = 9 (i“l,2, . . . ,n)  (4) 

and 

^(A-Xjl+BGC)  - 0 <j-l,2,...,n)  (5) 

so  that 


[A-Xil  , B] 


Si 

t“iJ 


and 


- s'] 


"--j 


where 


(i*l,2, . . . ,n) 


(6) 


(j-1,2,... ,n)  (7) 


“i 

(i«l,2, . . . ,n) 

(8) 

5j 

- G'B'$. 

- - - J 

( j“l » 2 , • • , ,n) 

(9) 

and 

" 6ij 

( i » J“1 1 2 , . . . ,n) 

. (10) 
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Conversely,  if  equations  (6),  (7),  and  (10)  are  satisfied  by 
a self-conjugate  set  {X^, X2 » • . . , X } of  distinct  complex 
numbers  and  corresponding  self-conjugate  sets  (x^, X2 * • • • • Xn) 
and  * * * '$n^  of  linearlY  independent  vectors,  then 

equations  (8)  and  (9)  are  satisfied  by  a matrix  G £ 
such  that  Xj , • • • » XR}  is  the  eigenvalue  spectrum  of  the 
matrix  (A+BGC)  and  { x, , X, » • • • » X„ } and  } constitute 

corresponding  eigenvector  and  reciprocal  eigenvector  sets. 

It  accordingly  follows  from  equations  (8)  and  (9),  respectively, 
that  the  real  output-feedback  matrix 

G « [u^  , w2  , , Up]  [CXj,  , CX2  , ...  , CXp] 

(11) 

and  the  real  transposed  output- feedback  matrix 

* ^5l  ' 5r2  ' ’ * ’ ' C?f$i  • ' ***  ' ^^nyl 

(12) 

assign  the  self-conjugate  distinct  eigenvalue  spectrum 
{X^, X2 , . . . , Xfi}  and  the  corresponding  eigenvector  and  reciprocal 
eigenvector  sets  (Xj/  X2 / • • . , Xn>  and  ^ t0  the 

closed-loop  plant  matrix  (A+BGC)  in  case  {Cx, ,CX^ , . . . ,CX  } 

~ ~ I — ~ 2 ~ ~ p 

is  a set  of  p linearly  independent  vectors  and  {B'$.,B'<J>  , 

— •►I  *•  ^2 

...,B'$m}  is  a set  of  m linearly  independent  vectors.  Such 
seta  (CXj^CXj, . . . ,CXp}  and  $2 ,... ,B* clearlY 

exist  when  rank  C « p,  rank  B * m,  and  {X,,X-,...,X  } and 
* * * r-n^  ar®  sets  linearly  independent  vectors. 

It  is  thus  evident  that,  in  the  case  of  self-conjugate 
distinct  eigenvalue  spectra,  closed-loop  eigenstructure  is 
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assignable  by  output  feedback  just  in  case  the  eigenvalue 
spectrum  {X^, *2* • • • * XR}  is  such  that  the  corresponding  eigen- 
vector and  reciprocal  eigenvector  sets  {X, , X. , . . . , X } and 

- 1 ~ i - n 

(♦  i»$2'  * * ’ '~n^  the  su^sPaces  determined  (Porter  and 

D'Azzo  1977)  in  accordance  with  equations  (6)  and  (7),  respect- 
ively, by  the  kernels  of  each  of  the  n matrices 


S(Xi)  - [A-Xjl  , B]  (i*l, 2 , . . . ,n)  (13) 

together  with  the  kernels  of  each  of  the  n matrices 

T'(Xj)  - [A'-Xjl  , C'J  (j«l,2,...,n)  . (14) 

3.  Illustrative  Example 

These  results  can  be  conveniently  illustrated  by 
designing  an  output-feedback  regulator  for  the  system 
governed  by  the  respective  state  and  output  equations 
(Davison  and  Wang  1975) 


x(t)  * 

'0,1,0' 

0,0,1 

x ( t ) + 

'1  , O' 

1 , 0 

0 

<* 

0 

-1  ' \ 

u(t) 


(15) 


and 


*(t) 


1,0,0' 
[0  , l , 0 


x(t) 


(16) 


such  that  the  eigenvalue  spectrum  of  the  closed-loop  plant 
matrix  is 


(J  ( A+BGC) 


I ^i , X2 , X3 } m {-1,-2,— 5} 


• (17) 


Indeed,  it  is  evident  from  equations  (13)  , (14) , (15) , and 
(16)  that 
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“A  f 1 f OflfO 

0 i ~ A f 1 i 1 i 0 

.0  , 0 , -X  , 1 , 1 


“A,  0 , 0 , 1 , Ol 


T'(A)  - 1 f -X  , 0 , 0 , 1 

.0  f X , -X  , 0 , O 


It  therefore  follows  immediately  (Porter  and  D'A22o  1977) 
from  equation  (18)  that 


ker  S(-l) 


span 


ker  S (-2) 


span 


ker  S (-5) 


span 


and  similarly  (Porter  and  D'Azzo  1977)  from  equation  (19) 
that 


ker  T' (-1) 


span 


ker  T' (-2)  « span 


ker  T' (-5)  * span 


It  is  thus  evident  from  equations  (20) , (21) , and  (22)  that 
the  closed-loop  eigenvectors  corresponding  to  the  eigenvalue 
spectrum  (17)  must  be  assigned  to  the  respective  subspaces 


26 
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Z(- 2) 


span 


E(-5) 


span 


and  it  is  similarly  evident  from  equations  (23) , (24)  , and 
(25)  that  the  closed-loop  reciprocal  eigenvectors  correspond- 
ing to  the  eigenvalue  spectrum  (17)  must  be  assigned  to  the 
respective  subspaces 


n-u 


span  0 , 1 


r (-2) 


span  0 , -2 


T(-5) 


span  0 , -5 


Since  the  vectors 
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-1  £ E(-1) 


2 E Z(-2) 


1 E ^(-5) 


3/4  £r(-i) 

-3/4 


2/3  E n-2) 


$3  = 5/12  En-5) 


1/12 


are  clearly  such  that  {x1#X2/X3>  and  {♦1#$2#$3}  constitute 
sets  of  linearly  independent  vectors  with  the  property  that 


*j*i 


(i,j=l,2,3)  , (38) 


it  follows  from  equations  (11)  and  (12)  that  equation  (17) 
is  satisfied  by  the  output-feedback  matrix 


The  corresponding  output-feedback  regulator  is  accordingly 
governed  by  the  control-law  equation 


-4  , -4 

u(t)  = y ( t)  . (40) 

-10  , -9J ' 

4.  Conclusion 

In  this  paper,  the  method  of  entire  eigenstructure  assign- 
ment has  been  applied  to  the  design  of  linear  multivariable 
continuous-time  output-feedback  regulators.  It  has  been  shown 
that,  in  the  case  of  self-conjugate  distinct  eigenvalue  spectra, 
the  closed-loop  eigenstructure  assignable  by  output  feedback  is 
constrained  by  the  requirement  that  the  elements  of  the  sets 
of  linearly  independent  self-conjugate  vectors  {X.^ , X2 » • • • » Xn) 
and  {«^,  $2'  * * * '$n^  lie  in  subsPaces  determined  by  the  kernels 
of  S(Xi)  (i=*l,2, . . . ,n)  and  T'(Xj)  ( j = l , 2 , . . . ,n)  , respectively, 
and  satisfy  the  orthogonality  conditions  (10) . In  constrast, 
the  closed-loop  eigenstructure  assignable  by  state  feedback  is 
constrained  only  by  the  requirement  that  the  elements  of  the 
set  of  linearly  independent  self-conjugate  vectors  { , X2 » • • • , Xn } 
lie  in  subspaces  determined  by  the  kernels  of  S ( X±)  (i=l,2, . . . ,n) . 

It  is  because  of  the  severe  constraints  on  the  closed-loop 
eigenstructure  assignable  by  output  feedback  that  it  is  frequently 
impossible  to  achieve  satisfactory  closed-loop  behaviour  by 
means  of  static  continuous-time  output-feedback  regulators,  and 
that  it  is  consequently  necessary  to  introduce  dynamic  compensator; 
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(Brasch  and  Pearson  1970) . However,  the  design  of  such  dynamic 
continuous-time  output-feedback  regulators  can  be  effected  by 
applying  the  method  of  entire  eigenstructure  assignment  in  the 
manner  of  Section  2 to  appropriately  augmented  (Brasch  and 
Pearson  1970,  Kimura  1975)  continuous-time  systems. 
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Abstract 

In  this  paper,  the  method  of  entire  eigonstructuro 
assignment  (Porter  and  D'Azao  1978a, b)  is  applied  to  the 
design  of  linear  multivariabLe  discrete-time  output-feedback 
regulators.  It  is  shown  that,  in  the  case  of  self-conjugate 
eigenvalue  spectra,  the  closed-loop  eigenstructure  assign- 
able by  output  feedback  is  constrained  by  the  requirement 
that  the  eigenvectors  and  generalised  eigenvectors  and  the 
reciprocal  eigenvectors  and  generalised  reciprocal  eigen- 
vectors lie  in  well-defined  subspaces.  The  method  is  illu- 
strated by  designing  an  output- feedback  regulator  for  a thii'd- 
order  discrete- time  system. 
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1.  Introduction 

In  this  paper,  the  method  of  entire  eigenstructure 
assignment  (Porter  and  D’Azzo  1978a, b)  is  applied  to  the 
design  of  output-feedback  regulators  for  multivariable  linear 
discrete- time  systems  governed  by  state  and  output  equations 
of  the  respective  forms 

x(k+l)  - Ax(k)  + Du (k ) (1) 

and 

y(k)  - Cx  (k)  , (2) 

where  a£  Rnxn,  B ££  Rnxm,  C ££  RpXn,  rank  B = m,  and  rank 
C * p.  The  theory  is  analogous  to  that  developed  by  Porter 
and  Bradshaw  (1978)  for  continuous-time  regulators  but  is 
significantly  extended  in  order  to  allow  the  assignment  of 
confluent  eigenvalues  to  the  plant  matrix  of  the  closed-loop 
system.  It  is  therefore  possible,  for  example,  to  apply  this 
theory  to  the  design  of  output- feedback  regulators  with 
finite  settling  times. 

2,  Theory 

Thus,  if  output  feedback  is  applied  to  the  system  governed 
by  the  state  equation  (1)  in  accordance  with  the  control-law 
equation 

u(k)  <*  G^(k)  (3) 

and  the  output-feedback  matrix  G £ Rmxp  is  such  that  the 
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closed-loop  plant  matrix  (A+BGC)  has  a self-conjugate  eigen 
value  spectrum  {A^, X2 , . . . , Afc} , a corresponding  eigenvector 
and  generalised  eigenvector  set(xx  :h»l,2, . . . ,111^; 

j«l,  2, . . . . . . , t} , and  a corresponding  reciprocal 

eigenvector  and  generalised  reciprocal  eigenvector  set 
{iji^  * i csl/  2 / • * * / mj3a;  t)“  1 / 2 / • • • a=  1 f 2 t • • • 1 1 } , then 

cl 

(A-X.I+BGC)  (1'3)  - 0 (4a) 

(A-X.I+BGC) x, (2' j)  » Xl<l'3>  (4b) 

~ 1 Ai  ~Ai 


(m..-l,j) 

<*-*!**«>  X*i  = XXi31 

(jssl»2,...,k.;  i®!  , 2 , . . . , t)  t 


(4m;ji) 


and 


(1'b)  ' (A-X.I+BGC)  - 4)  (2'b)' 


a~  ~ ~ ~ 


XX 


(5a) 


♦x 


(mba~1,b)  ’ 


(A-X.I+BGC)  = <p 


(mhatbr 


— «•  * 


z\ 


a 


(mba,b> 

<J>,  Da  (A-X  I+BGC)  - 0' 

^ A _ cl  -v  ^ ^ «• 


(5"W 


( h= 1,2, • • • /k^; a* 1 , 2 , • • • , t)  , 

where  X\^'^  is  the  hth  vector  in  the  jth  string  of  length 
~ Ai 

associated  with  the  eigenvalue  Xi#  and  where  $^c'b^  is 
the  cth  vector  in  the  bth  string  of  length  m^a  associated 
with  the  eigenvalue  A . The  vectors  ( j-1,2, . . . ,k. ) 

a A ^ 1 
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b 


i 


are  the  k^  eigenvectors  associated  with  the  eigenvalue 

whilst  the  remaining  vectors  in  each  of  the  ki  strings  of 

vectors  satisfying  equations  (4)  are  generalised  eigenvectors 

associated  with  the  eigenvalue  A. . Similarly,  the  vectors 
(mba'b) 

<j>,  (b=l,2, . . . ,k  ) are  the  k reciprocal  eigenvectors 

a a a 

a 

associated  with  the  eigenvalue  A , whilst  the  remaining 

a 

vectors  in  each  of  the  k strings  of  vectors  satisfying 

a 

equations  (5)  are  reciprocal  generalised  eigenvectors  associated 

with  the  eigenvalue  A . The  total  number  of  vectors  associated 

a 

with  the  eigenvalue  A^  in  each  set  is  evidently 


mf 


Z m 
g=i 


gf 


(f=l,2,...,t) 


(6) 


and 


t 

n = Em-  (7) 

f=l  1 


Equations  (4)  and  (5)  can  be  written  in  the  form 


&-X±l 


[A-AiJ 


?] 


£a  . 


(1/j) 

i 


a 


u. 


5A 


S?A 


dr  j) 
'i 

(2,  j) 
i 

(2 , j) 


0 


(8a) 


9 


(8b) 
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(c,b) ' (h,j) 

Z\. 

a i 


6ai5bj6ch 


(h=l / 2 , < . • , m j ^ ; j*l , 2 , . . . ,k^ji— 1,2, • . . , t) 

(c-1,2,... ,mba;b=L,2,...,ka;a=l,2,.q.,t)  . 

Conversely,  if  equations  (8),  (9),  and  (12)  are  satisfied 
by  a self-conjugate  set  { , Xj , . • • , Afc}  of  complex  numbers 
and  corresponding  self-conjugate  «5ets  { X,  :h»l , 2 , . . . ,m  . . ; 

- J 1 

j *1  # 2 * • • • • • • * t } and  { r ^ i c53!  f 2 9 • • • t J 

a 

b=l,2, . . . ,ka;a=l,2 , . . . ,t}  of  linearly  independent  vectors, 
then  equations  (10)  and  (11)  are  satisfied  by  an  mxn  matrix 
G such  that  (A^, A^ » . . . , Afc}  is  the  eigenvalue  spectrum  of  the 


matrix  (A+BGC) , (x,  (h' ^ sk*l, 2 , . . . ,m.  . ? j-1,2 , . . . ,k . ; 

- ---  A ^ 31  1 

i*l,2,...,t}  constitutes  a corresponding  eigenvector  and 

generalised  eigenvector  set,  and  { (f ^ : c=l , 2 , . . . ,mba; 

a a 

b“l,2, . . . ,ka;a=l,2, . . . ,t}  constitutes  a corresponding  reciprocal 
eigenvector  and  reciprocal  generalised  eigenvector  set.  It 
accordingly  follows  from  equations  (10)  and  (11)  respectively 
that  the  real  output- feedback  matrix 


:k*l,2, . . . , m j i?j=*l»2,.  • • » k i ; 


QSi  * 


$ * * • r U> 


p^  f-9^1  ' S*2  ' * " ’ ' ?^p] 


and  the  real  transposed  output- feedback  matrix 


fii  ' 5; 


sJS'h  • r*2 


B'$  T 
- imJ 


(14) 


Indeed,  it  is  evident  from  equations  (13),  (14),  (15),  and 
(16)  that 


assign  the  self-conjugate  eigenvalue  spectrum  {X^, X^, . . . , Xfc}, 
the  corresponding  eigenvector  and  generalised  eigenvector  set 
{ ^ :hml, 2, . . . , m j ^jj—l,2, • • • ,k^; i = 1 , 2 , • • • ,t) , and  the 
corresponding  reciprocal  eigenvector  and  reciprocal  generalised 
eigenvector  set  (C,b* !C-1, 2, . . . ,mba}b=l, 2, . . . ,ka?a»l, 2, . . . , t} 

cl 

to  the  closed-loop  plant  matrix  (A+BGC)  in  case  (CX, ,CX- , . . . ,CX  } 

~~1  --2  ~-p 

is  a subset  of  p linearly  independent  members  of  the  set  (CX.^»j). 

~-Ai 

h**l ,2,  * • • , m j ^ j j“l, 2 , • • • ,k^ ; i=l , 2 , « • • , t } and  • • • , B * 

is  a subset  of  m linearly  independent  members  of  the  set 
!C-1,2, . . . ,mba;b=l,2, . . . , ka<*  a-1, 2 , . . . , t} . 

a 

It  is  thus  evident  that  closed-loop  eigenstructure  is 
assignable  by  output  feedback  just  in  case  the  self-conjugate 
eigenvalue  spectrum  (X^, X2, . . . , Xfc)  is  such  that  the  correspond- 
ing eigenvector/generalised  eigenvector  and  reciprocal  eigen- 
vector/reciprocal generalised  eigenvector  sets  {X^h'^: 
h* 1# 2, * ■ • , m j ^ ; j=l, 2, . . . ,k^ ; i® l,2,**»,t}  and  { $ ^ 9 * c* 1,2 , ■ • t , m^a ^ 
b*l, 2 , . . . ,ka; a*l, 2 , . . . , t)  lie  in  the  subspaces  determined 
(Porter  and  D'Azzo  1978b)  in  accordance  with  equations  (8)  and 
(9) , respectively,  by  the  kernels  and  generalised  kernels  of 
each  of  the  t matrices 


{5'h(c'b)" 


S(Xi) 


[A-Xjl  , B] 


( i=l ,2, • • • ,t) 


(15) 


together  with  the  kernels  and  generalised  kernels  of  each  of 
the  t matrices 

T'(Xa)  - [A'-XaI  , C']  (a-l,2,...,t)  . (16) 

3,  Illustrative  Example 

These  results  can  be  conveniently  illustrated  by  designing 


an  output- feedback  regulator  for  the  system  governed  by  the 


respective  state  and  output  equations 


x(k+l) 


and 


*(k) 


0,1,0' 

0 , o' 

- 

1,1,0 

x(k)  + 

1 , 0 

0,0,1 

H 

o 

o 

s 

x (k) 

0,0,  1 

u(k) 


(17) 


(18) 


such  that  x(k)  =0  after  a finite  number  of  discrete-time 
intervals.  Indeed,  it  is  evident  from  equations  (15),  (16), 
(17) , and  (18)  that 


S ( X) 


" A , 1 , 0 ,0,0 

1 , 1-X  , 0 ,1,0 

0 , 0 , 1-X  , 0 , 1 


and 


(19) 


(20) 


In  this  case,  it  is  necessary  to  assign  the  value  zero  to 
all  the  eigenvalues  of  the  closed-loop  plant  matrix.  Such 
an  assignment  clearly  requires  that  t - 1,  - 3,  and  X^  - 0 

and  therefore,  in  consonance  with  the  results  of  Rosenbrock 
and  Hayton  (1977),  that  k « 1 and  mn  = 3.  It  therefore 
follows  immediately  (Porter  and  D'Azzo  1978b)  from  equation 
(19)  that 


1 


I 1 

1 ( 


and  similarly  (Porter  and  D'Azzo  1978b)  from  equation  (20) 
that 


(22) 


It  is  evident  from  equation  (21)  that  the  closed-loop  eigen- 
vector corresponding  to  the  eigenvalue  »=  0 must  be 
assigned  to  the  subspace 


whilst  the  associated  string  of  generalised  eigenvectors 

Xq  1 ^ and  X^3,1)  must  be  generated  in  accordance  with  equations 

(8) , and  it  is  similarly  evident  from  equation  (22)  that 

the  closed-loop  reciprocal  eigenvector  corresponding  to 

the  eigenvalue  X^  = 0 must  be  assigned  to  the  subspace 


u 

0 
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span  < -1  , 0 


whilst  the  associated  string  of  generalised  reciprocal  eigen- 
vectors ^ and  roust  be  generated  in  accordance  with 


equations  (9) . Since  the  vectors 


yU,D 

-o 

J1'1* 

-o 


Y (2  # 1) 
io 


(3,1) 


J3'1) 

-O 


-o 


- (3,1) 
-o 


o*2'1* 

Zo 


(2,1) 


and 


1 

0 

-1 

-1 

2 


‘o' 

to 

0 

a 

1 

,(1,1) 

_-o 

1 

_-2_ 
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are  clearly  such  that 


, (c,l) \(h,l) 
?o  *o 


(29) 


(30) 


ch 


(c,h=l, 2, 3) 


(31) 


it  follows  from  equations  (13)  and  (14)  that  the  required 
eigenstructure  assignment  is  achieved  by  the  output  feedback 
matrix 


-2  , 1* 
1 , "2 


(32) 


The  corresponding  output-feedback  regulator  is  accordingly 
governed  by  the  control-law  equation 


u(k) 


-2  , 1 

1 , -2 


y (k) 


(33) 


It  can  be  readily  verified  that  the  state  of  the  closed-loop 
system  governed  by  equations  (17),  (18),  and  (33)  is  reduced 
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from  any  initial  value  to  zero  in  at  most  three  discrete- 
time intervals,  as  required. 


4,  Conclusion 


In  this  paper,  the  method  of  entire  eigenstructure 
assignment  has  been  applied  to  the  design  of  linear  multi- 
variable  discrete-time  output  feedback  regulators.  It  has 
been  shown  that  the  closed-loop  eigenstructure  assignable  by 
output  feedback  is  constrained  by  the  requirement  that  the 
elements  of  the  sets  of  linearly  independent  self-conjugate 
vectors  {X  :h=l,2, . . . ,m. . ; j=l, 2 , . . . ,k . ; i*l,2, . . . , t } and 

Jl  1 

!C=1,2»  • • • ,n*ba;b=l,2, . . . ,ka;a=l,2, . . . ,t}  lie  in  suo- 
a 

spaces  determined  by  the  kernels  and  generalised  kernels  of 
S(L)  (i=l,2, . . . ,t)  and  T'(X  ) (a=»l,2, . . . ,t)  , respectively,  and 
satisfy  the  orthogonality  conditions  (12) . In  contrast,  the 
closed-loop  eigenstructure  assignable  by  state  feedback  is 
constrained  only  by  the  requirement  that  the  elements  of  the 
set  of  linearly  independent  self-conjugate  vectors 

~Ai 

h*l,2, . . . ,m^j  j»l,2, . . . ,k^;i*l,2, . . . ,t)  lie  in  subspaces 
determined  by  the  kernels  and  generalised  kernels  of  S(X^) 
(i=l,2, . . . , t) . It  is  because  of  the  severe  constraints  on 
the  clor.ed-loop  eigenstructure  assignable  by  output  feedback 
that  it  is  frequently  impossible  to  achieve  satisfactory 
closed-loop  behaviour  by  means  of  static  discrete-time  output- 
feedback  regulators,  and  that  it  is  consequently  necessary  to 
introduce  dynamic  compensators  (Brasch  and  Pearson  1970) . How- 
ever, the  design  of  such  dynamic  discrete-time  output-feedback 
regulators  can  be  effected  by  applying  the  method  of  entire 


eigenstructure  assignment  in  the  manner  of  Section  2 to 
appropriately  augmented  (Brasch  and  Pearson  1970)  discrete- 
time systems. 
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APPENDIX  7 


DESIGN  OF  LINEAR  MULTIVARIABLE 
CONTINUOUS-TIME  TRACKING  SYSTEMS 
INCORPORATING  ERROR-ACTUATED  DYNAMIC  CONTROLLERS 


Professor  B Porter  and  Dr  A Bradshaw 

Department  of  Aeronautical  and  Mechanical  Engineering 
University  of  Salford 
Salford  M5  AWT 
England 


(International  Journal  of  Systems  Science,  vol  9,  pp627-637,  1978) 


Abstract 

In  this  paper,  the  method  of  entire  eigenstructure 
assignment  (Porter  and  Bradshaw  1978)  is  applied  to  the 
design  of  linear  multivariable  continuous-time  tracking 
systems  incorporating  error-actuated  dynamic  controllers. 

The  method  is  illustrated  by  designing  an  error-actuated 
dynamic  controller  which  causes  the  output  of  a second- 
order  continuous-time  plant  to  track  a constant  command 
input  in  the  presence  of  an  unmeasurable  constant  disturbance 
input. 
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1.  Introduction 


In  this  paper,  the  method  of  entire  eigenstructure 
assignment  (Porter  and  Bradshaw  1978)  is  applied  to  the 
design  of  linear  multivariable  continuous -time  tracking 
systems  incorporating  error-actuated  dynamic  controllers. 
Such  tracking  systems  consist  of  a controllable  and  observ- 
able nth-order  linear  multivariable  plant  governed  by  state 
and  output  equations  of  the  respective  forms 


x(t)  = Ax(t)  + Bu  ( t)  + Dd(t) 


y ( t)  = Cx(  t) 


where  B and  C have  full  rank,  together  with  an  error-actuated 
dynamic  controller  which  is  required  to  cause  the  pxl  output 
vector  y(t)  of  the  plant  to  track  a pxl  command  input  vector 
r(t)  in  the  sense  that 


lim  e(t)  = lim{r (t) -y (t) } = 0 

t-**®  ~ t-*00  ~ 

for  polynomial  command  and  disturbance  inputs  of  the 
respective  forms 


r ( t) 


I 6,  ,tJ 
i=l  1 1 


. IH 
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It  is  important  to  note  that  tracking  systems  incorporating 
error-actuated  dynamic  controllers  can  be  designed  for  a 
much  larger  class  of  plants  than  tracking  systems  incorporat- 
ing error-actuated  static  controllers  (Porter  and  Bradshaw 
1976)  in  view  of  the  fact  that  eigenstructure  assignment  by 
error-actuated  static  controllers  and  by  output-feedback 
controllers  (Kimura  1975,1977)  are  essentially  equivalent. 

2.  Theory 

The  first  stage  in  the  design  of  the  required  error- 
actuated  dynamic  controller  for  the  plant  governed  by 
equations  (1)  and  (2)  involves  the  introduction  (Porter  and 
Bradshaw  1974)  of  a vector  comparator  and  a series  of 
q = max(r,s)  vector  integrators  in  order  to  generate  the  q 
vectors  defined  by  the  equations 

^(t)  * e ( t)  , 

z2(t)  - r 1 ( t ) , 


?q(t)  " 5q-l(t)  * j 

It  is  then  evident  from  equations  (1),  (2),  and  (5)  that 
the  open-loop  tracking  system  is  governed  by  state  and 
output  equations  of  the  respective  forms 


and 


The  open-loop  tracking  system  governed  by  equations  (7) 
and  (8)  is  controllable  if  and  only  if  (Porter  and  Bradshaw 
1974) 


rank 


n + p 


(9) 


since  (A,B)  and  (£,&)  are  respectively  controllable  and 
observable  pairs. 
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In  the  case  of  such  controllable  and  observable  open- 
loop  tracking  systems,  the  second  stage  in  the  design  of  the 
error-actuated  dynamic  controller  involves  the  introduction 
of  an  ith-order  dynamic  compensator  (Brasch  and  Pearson  1970) 
governed  by  state  and  output  equations  of  the  respective  forms 


w(t) 


Fw(t)  + Ge ( t)  + E H.z. (t) 

l-l'i-l 


Kw(t)  + Le (t ) + E M.z. (t) 


where 


min(vc-l,vo-l) 


and  vc  and  are  respectively  the  controllability  and 
observability  indices  of  the  open-loop  tracking  system 
governed  by  equations  (7)  and  (8) . It  is  then  evident  from 
equations  (7),  (8),  (10),  and  (11)  that  the  closed-loop 
tracking  system  is  governed  by  state  and  output  equations 
of  the  respective  forms 


*l(t) 


!q(t) 


"A-BLC  , BMX 
-C  , 0 


0 , 0 
-GC  , H, 


BM  , BKl  Tx(t)  ' 
0 , 0 z2(t) 


0 , 0 zq(t) 


Uq  * !, 


w(t) 
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+ r (t)  + d(t) 


y(t) 

5i(t) 


!q<t) 

w(t) 


C , 0 , 

9 ' Jp  ' 


0,0, 

0,0, 


0 , 0 ’x(t) 

0 , 0 z-^t) 


• Jp  t 9 

, 0 , IA  w(t) 


. (14) 


In  view  of  equations  (4)  and  (5) , it  is  clear  by  differentiat- 
ing equation  (13)  (q-1)  times  that  the  closed-loop  tracking 
system  will  behave  so  that 


= 0 


and  therefore,  in  view  of  equations  (6) , so  that  equation 
(3)  will  be  satisfied  if  the  error-actuated  dynamic  com- 
pensator governed  by  equations  (10)  and  (11)  is  designed 
such  that  all  the  eigenvalues  of  the  plant  matrix  of  the 
closed-loop  tracking  system  governed  by  equations  (13)  and 
(14)  are  assigned  to  the  open  left-half  of  the  complex  plane. 


3.  Compensator  Design 


It  is  evident  from  equations  (7),  (8),  (10),  and  (11) 
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that  such  a compensator  can  be  designed  by  the  synthesis  of 
an  appropriate  output-feedback  control  law  of  the  form 


Thus,  if  the  (m+£) x(p+pq+£)  output- feedback  matrix  G.  is 

% X 

such  that  the  closed-loop  plant  matrix  (A^+B^G^C^)  has  a 
self-conjugate  spectrum  (X1# X2, • • • * xn+pq+fc*  of  distinct 
eigenvalues  and  corresponding  eigenvector  and  reciprocal 
eigenvector  sets  . . . <Xn+pq+4}  and  • • • »$n+pq+ih 

then  obviously 

(Ajl-XiI+BjlGJlC4)xi  * 0 (i-l,2,...,n+pq+i)  (26) 


ai. 


and 
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« 0 

•m 

(j-1,2,. . . , n+pq+l) 

(27) 

so  that 

0 (i«l,2, . . . , n+pq+l)  (28) 


0 (j-l»2f... , n+pq+l)  (29) 


where 


“i  = 

SiStXi 

(i-1,  2 , . . . , n+pq+l) 

(30) 

5j  3 

Sl?i*  j 

(j-1,2, . . . , n+pq+l) 

(31) 

and 

♦5*i 

* 6ij 

(i,  j-1,2, . . ., n+pq+l) 

. (32) 

Conversely,  if  equations  (28) , (29)  , and  (32)  are  satisfied 
by  a self-conjugate  set  (X^, X2, . . . , *n+pq+j^  of  distinct 
complex  members  and  corresponding  self-conjugate  sets 

{*1'*2 *nWt>  and  (Jj.+j *nwl)  °f  linearly 

independent  vectors,  then  equations  (30)  and  (31)  are  satisfied 
by  an  (m+1)  x (n+pq+l)  matrix  G£  such  that  { X 1 , X2 , • • • » *n+pq+e * 
is  the  eigenvalue  spectrum  of  the  matrix  (A.+B.G.C. ) and 
{XiaX2f*»Xn+pq+Jl>  and  U1'!2'**-'$n+pq+£}  constitute 
corresponding  eigenvector  and  reciprocal  eigenvector  sets,  it 
accordingly  follows  from  equations  (30)  and  (31)  respectively 
that  the  real  output- feedback  matrix 
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-8.  “ fei  ' ~2  ' * * ’ ' 9p+pq+*3  ' 9**2  ' * * * ' 9i*p+pq+^ 


-1 


(33) 


and  the  real  transposed  output- feedback  matrix 


s; 


fei  ' 5: 


9 • • • 9 


5m+J  ftih  - 


• » B'J 


-1 


riirm+ji1 


(34) 


assign  the  self-conjugate  distinct  eigenvalue  spectrum 

(X^  *2'  * * * '*n+pq+Jt^  and  the  corresPondin9  eigenvector  and 
reciprocal  eigenvector  sets  { Xx , X2 • • • • * Xn+pq+i ) and 

^$1'$2' * * * '^n+pq+4^  to  the  closed-looP  plant  matx  u:  (A^+B^G^C^) 
in  case  t£j!,£l' 9)t?2  ' • * • '9s3p+pq+j^  is  a set  of  (p+pq+£) 
linearly  independent  vectors  and  (B'4^  ,B'$2 , . . . , B^  A } is 
a set  of  (m+i)  linearly  independent  vectors,  respectively. 

In  view  of  equations  (28),  (29),  (33),  and  (34),  the 
computation  of  is  thus  reduced  to  the  determination 
(Porter  and  D'Azzo  1977)  of  the  kernels  of  each  of  the  n 


matrices 

S£  ( A±)  = [A^-X.^1  , Bj  ( i=l , 2 , . . . , n+pq+X, ) (35) 

together  with  the  kernels  of  each  of  the  n matrices 

= [A^X.I  , C (j-1,2,..., n+pq  + 4.)  (36) 

followed  by  the  selection  of  sets  of  linearly  independent 


self-conjugate  vectors  (x^Xj Xn+pq+ll>  and  *n+pq+t> 

from  subspaces  determined  by  the  kernels  of  ( A ^) (i»l, 2, . . . ,n+pq+£; 
and  T' ( X j ) (j=l,2,... , n+pq+fc) , respectively,  such  that  the 


pFTTT— m — m — IMS , 
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orthogonality  conditions  (32)  are  satisfied.  It  is  finally 
evident  from  equations  (10),  (11),  and  (22)  that  the  matrices 
in  the  respective  state  and  output  equations  of  the  required 
£th-order  dynamic  compensator  are  determined  by  the  sub-matrices 
of  the  output-feedback  matrix  G^. 
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4.  Illustrative  Example 


The  results  presented  in  Sections  2 and  3 can  be 
conveniently  illustrated  by  designing  an  error-actuated 
dynamic  controller  which  will  cause  the  output  of  the 
controllable  and  observable  linear  plant  governed  by  the 
respective  state  and  output  equations 


x^t) 
x2  (t) 


o , ii  rXl(t) 


i , ij  Lx2(t) 


+ 1 u(t)  + I d (t)  (37) 


y (t) 


x,  (t) 

[1.0]  1 


to  track  any  constant  command  input  r(t)  in  the  presence 
of  any  unmeasurable  constant  disturbance  input  d(t).  In 
this  case  it  is  clear  that  r *=  s = q = 1,  so  that  the 
open-loop  tracking  system  is  governed  by  the  respective 
state  and  output  equations 


xi(t) 
x2  (t) 


0 , i , oi  rx1(t)-|  (o' 

1 • 1 , 0 x2(t)  + 1 u(t) 

-1,0,0  2 ( t ) 0 


+ 0 r (t)  + Md(t) 


i 
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y (t) 


1,0/0 

.0,0,1 


Xj^t) 

x2(t) 


This  system  is  clearly  controllable  and  observable  with 
\>c  - 3 and  vQ  = 2 so  that  (Brasch  and  Pearson  1970)  equation 
(12)  indicates  that  i = 1.  Furthermore,  in  the  notation  of 
equations  (23),  (24),  and  (25),  it  follows  from  equations 
(39)  and  (40)  that 


0 ,1,0,0 

1 ,1,0,0 

-1  , 0 , 0 , 0 
0 ,0,0,0 


0 , 0 

1 , 0 

o , 0 

! 0 , 1 


1 , 0 , 0 , 0 

0 , 0 , 1 , 0 

0 , 0 , 0 , 1 


It  is  thus  evident  from  equations  (35),  (36),  (41),  (42), 
and  (43)  that 


?llX) 


"X  , 1 ,0, 0,0,0 

1 , 1-X  , 0 , 0,1,0 

**  1 , 0 , -X  , 0,0,0 

o , 0 , 0 , “*  X ,0,1 


In  order  to  design  an  error-actuated  dynamic  compensator 
for  th'i  open- loop  tracking  system  governed  by  equations 
(39)  and  (40)  such  that  the  eigenvalues  of  the  plant  matrix 
of  the  closed-loop  tracking  system  are 


UlfX2,X3,X4)  = (-1,-2, -3, -4}  (46) 

the  design  method  described  in  Section  3 can  be  readily 
used  to  compute  an  output-feedback  matrix  such  that 

o(£l+?l9l9l)  = (-1,-2, -3, -4}  . (47) 


Indeed,  it  follows  immediately  (Porter  and  D’Azzo  1977) 
from  equation  (44)  that 


(48) 


AD-A071  662 


UNCLASSIFIED 
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It  is  evident  from  equations  (48),  (49),  (50),  and  (51)  that 
the  closed-loop  eigenvectors  corresponding  to  the  eigenvalue 
spectrum  (47)  must  be  assigned  to  the  respective  subspaces 


and 


I 
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and  it  is  similarly  evident  from  equations  (52),  (53),  (54), 
and  (55)  that  the  closed-loop  reciprocal  eigenvectors 


corresponding  to  the  eigenvalue  spectrum  (47)  must  be 


?3 


and 


U 
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-16 

4 

13 

5 


€ ri<“3> 


35/6 

-7/6 

■14/3 

-V3 


€ V-4) 


art  clearly  such  that 


i'ih  - {1J 


1 ’ 


(70) 


(71) 


(i, j-1,2,3,4) 


(72) 


it  follows  from  aquations  (33)  and  (34)  that  aquation  (47) 
is  satisfiad  by  tha  output- feedback  matrix 


-47  , 34  , 10 

G,  - . (73) 

[49  , -35  , -11 

In  viaw  of  aquations  (10),  (11),  (22),  (73),  tha  corraspond 
inq  dynamic  compansator  for  tha  opan-lcop  tracking  system 
govarned  by  aquations  (39)  and  (40)  is  governed  by  tha 
raspactiva  stata  and  output  aquations 

w(t)  - -llw(t)  - 49a (t)  - 35s (t)  (74) 

and 

u ( t ) - 10w(t)  ♦ 47a(t)  ♦ 34a (t)  , (75) 

so  that  tha  raquirad  arror-actuatad  dynamic  contrellar  is 
charactarisad  by  tha  transfar  function 
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T(s)  - u ( s ) /e ( s ) - (47s24-61s+24)/s(s+ll)  . (76) 

It  can  be  readily  verified  that  the  poles  of  the  closed- 
loop  tracking  system  governed  by  equations  (37),  (38), (74), 
and  (75)  are  (-1,-2, -3, -4}  and  that 

lim  e(t)  - lim{r(t) -y  (t) ) - 0 (77) 

t-K»  t-*-» 

for  any  constant  command  input  r(t)  and  any  constant  un- 
measurable disturbance  input  d(t) . 

5.  Conclusion 

In  this  paper,  the  method  of  entire  eigenstructure 
assignment  has  been  applied  to  the  design  of  linear  multi- 
variable  continuous -time  tracking  systems  incorporating 
error-actuated  dynamic  controllers.  It  has  been  indicated 
that  such  tracking  systems  can  be  designed  for  a much 
larqer  class  of  plants  than  tracking  systems  Incorporating 
error-actuated  static  controllers  (Porter  and  Bradshaw  1976) 
in  view  of  the  fact  that  eigenstructure  assignment  by  error- 
actuated  static  controllers  and  by  output-feedback  controllers 
(Kimura  1975,1977)  are  essentially  equivalent. 
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APPENDIX  8 


DESIGN  OF  LINEAR  MULTIVARIABLE 
DISCRETE-TIME  TRACKING  SYSTEMS 
INCORPORATING  ERROR-ACTUATED  DYNAMIC  CONTROLLERS 


Dr  A Bradshaw  and  Professor  B Porter 

Department  of  Aeronautical  and  Mechanical  Engineering 
University  of  Salford 
Salford  M5  4WT 
England 


(International  Journal  of  Systems  Science.,  vol  9,  pp  1079-1090,  1978) 


Abstract 


In  this  paper,  the  method  of  entire  eigenstructure 
assignment  (3radshaw  and  Porter  1978a)  is  applied  to  the 
design  of  linear  multivariable  discrete-time  tracking 
systems  incorporating  error-actuated  dynamic  controllers. 

The  method  is  illustrated  by  designing  an  error-actuated 
dynamic  controller  which  causes  the  output  of  a second- 
order  discrete-time  plant  to  track  a constant  command  input 
in  the  presence  of  an  unmeasurable  constant  disturbance  input. 


F! 
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1.  Introduction 

In  this  paper,  the  method  of  entire  eigenstructure 
assignment  (Bradshaw  and  Porter  1978a)  is  applied  to  the 
design  of  linear  multivariable  discrete-time  tracking  systems 
incorporating  error-actuated  dynamic  controllers.  Such 
tracking  systems  consist  of  a controllable  and  observable 
nth-order  linear  multivariable  plant  governed  by  state  and 
output  equations  of  the  respective  forms 

x(k+l)  ■ Ax(k)  + Bu(k)  + Dd(k)  (1) 

and 


y (k)  =■  Cx  (k)  (2) 

where  B and  C have  full  rank,  together  with  an  error-actuated 
dynamic  controller  which  is  required  to  cause  the  pxl  output 
vector  y(k)  of  the  plant  to  track  a pxl  command  input  vector 
r(k)  in  the  sense  that 
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The  theory  developed  in  this  paper  is  analogous  to  that 
developed  by  Porter  and  Bradshaw  (1978)  for  continuous-time 
tracking  systems  but  is  significantly  extended  in  order  to 
allow  the  assignment  of  confluent  eigenvalues  to  the  plant 
matrix  of  the  closed-loop  tracking  system.  It  is  therefore 
possible  to  apply  this  theory  to  the  design  of  an  error- 
actuated  dynamic  controller  which  causes  the  output  vector 
of  a plant  governed  by  equations  (1)  and  (2)  to  track  a 
command  input  vector  in  the  sense  that 

e(k)  = r(k)  - y(k)  = 0 (k=  v,  vfl , . . . ) (6) 

for  command  and  disturbance  inputs  defined  by  equations  (4) 
and  (5) , where  v is  the  index  of  nilpotency  of  the  closed- 
loop  plant  matrix  of  the  tracking  system.  It  is  important  to 
note  that  tracking  systems  incorporating  error-actuated 
dynamic  controllers  can  be  designed  for  a much  larger  class 
of  plants  than  tracking  systems  incorporating  error-actuated 
static  controllers  (Bradshaw  and  Porter  1978b)  in  view  of 
the  fact  that  eigenstructure  assignment  by  error-actuated 
static  controllers  and  by  output- feedback  controllers  (Kimura 
1975,1977)  are  essentially  equivalent. 

2.  Theory 

The  first  stage  in  the  design  of  the  required  error- 
actuated  dynamic  controller  for  the  plant  governed  by  equations 
(1)  and  (2)  involves  the  introduction  (Bradshaw  and  Porter 
1975)  of  a vector  comparator  and  a series  of  q = max(r,s) 
discrete-time  vector  integrators  in  order  to  generate  the  q 
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0 


y 


Li 


0 


vectors  defined  by  the  equations 


z1(k+l)  = z^k)  + e (k) 
z2(k+l)  = z2(k)  + zx(k) 


?q(k+1)  " + 5q-l(k} 


(7) 


It  is  then  evident  from  equations  (1),  (2),  and  (7)  that 
the  open-loop  tracking  system  is  governed  by  state  and 
output  equations  of  the  respective  forms 


x(k+l) 

(k+1) 

= 

5q-l(k+1) 

-?,<k) 

^ f Of...,  0,  0 

' Jp  ' • • • • 9 ' ® 


0 , 0 


> Jp  ■ 2 


9 • 


' b ■ “I 


x(k) 

zx(k) 

5q-l (k) 

and 


p — 

B 

'o' 

■p 

0 

I 

o 

«*» 

~p 

• 

• • 

• 

+ 

• 

u(k)  + 

• • 

r(k)  + 

• 

d (k) 

0 

0 

0 

*p 

•» 

0 

0 

o 

y(k) 

"c 

9 

0 , ... 

9 

0 , O' 

««P 

'x(k)  ' 

5i(k) 

0 

*>» 

9 

?p  • * * * 

9 

0 f 0 

p*p 

5i(k) 

= 

5q-l<k) 

0 

9 

0 f • • • 

9 

Jp  • ? 

5q-l(k> 

L*q<k>  J 

0 

9 

0 / • • • 

t 

9 ' Jp- 

-!q(k)  - 

(8) 


(9) 


LJ 
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The  open-loop  tracking  system  governed  by  equations  (8)  and 
(9)  is  controllable  if  and  only  if  (Bradshaw  and  Porter  1975) 


rank 


B 

■*» 

S 


A- 1 
-C 


n + p 


(10) 


since  ( A,  B)  and  (C,A)  are  respectively  controllable  and 
observable  pairs. 

In  the  case  of  such  controllable  and  observable  open- 
loop  tracking  systems,  the  second  stage  in  the  design  of  the 
error-actuated  dynamic  controller  involves  the  introduction 
of  an  ith-order  dynamic  compensator  (Brasch  and  Pearson  1970) 
governed  by  state  and  output  equations  of  the  respective 
forms 

q 

w(k+l)  - Fw(k)  + Ge(k)  + Z H . z,  (k)  (11) 

i-i'i'i 

and 

q 

u(k)  - Kw(k)  + Le  (k)  + Z M.  z. (k)  (12) 

~ --  i-l*'1'’1 

where 


l - min(vc-l,vo-l)  (13) 

and  vc  and  vQ  are  respectively  the  controllability  and 
observability  indices  of  the  open-loop  tracking  system 
governed  by  equations  (8)  and  (9).  It  is  then  evident  from 
equations  (8),  (9),  (11),  and  (12)  that  the  closed-loop 
tracking  system  is  governed  by  state  and  output  equations 
of  the  respective  forms 
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x(k+l) 

z,(k+l) 

**  X 


?q(k+l) 

w(k+l) 


A-BLC  , BM,  , 

“C  » In  * 


0,0, 

+0  «• 

“99  ' »i  * 


* BMq  , BK'  x(k) 
,0,0  5l(k) 


, I„  , 0 2 (k) 

~p  * -q 

, h_  , F.  w(k)  . 


r (k)  + d(k) 


y(k)  ' 
?l(k) 


z (k) 
~q 

w(k) 


C , 0 

“00 

0 ' 

~ -p 


0 , 0 

e*  •» 

0 , 0 


o , oi  rx(k) 

0 , O 2x(k) 


I_  , o z (k) 

-p  - ~q 


9 ' hi 


w(k) 


. (15) 


In  view  of  equations  (4)  and  (5),  it  is  clear  by  differencing 
equation  (14)  (q-1)  times  that  the  closed-loop  tracking  system 
will  behave  so  that 


lim  A ^ z (k)  ■ 0 

k—  “q 


and  therefore,  in  view  of  equations  (7),  so  that  equation  (3) 
will  be  satisfied  if  the  error-actuated  dynamic  compensator 
governed  by  equations  (11)  and  (12)  is  designed  such  that  all 
the  eigenvalues  of  the  plant  matrix  of  the  closed-loop  tracking 


J 
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system  governed  by  equations  (14)  and  (15)  are  assigned  to 
locations  within  the  unit  circle  of  the  complex  plane. 


3.  Compensator  Design 


It  is  evident  from  equations  (8) , (9) , (11) , and  (12) 
that  such  a compensator  can  be  designed  by  the  synthesis  of 
an  appropriate  output-feedback  control  law  of  the  form 

i?A(k)  * (17) 

for  the  augmented  open-loop  system  governed  by  state  and 
output  equations  of  the  respective  forms 


(k+1)  « AgX^tk)  + B^u^tk) 

and 


3ft  (k>  ' SiSt(k)  ' 


where 


i?*<k> 


u(k) 

v(k) 


?JL(k) 


x(k)  I 

!i(k> 


z (k) 

-q 

_w(k) 


(18) 


(19) 


(20) 


(21) 


Thus,  if  the  (m+i) x (p+pq+£)  output-feedback  matrix  G£  is 
such  that  the  closed-loop  plant  matrix  (A.+B-G.C, ) has  a 
self-conjugate  eigenvalue  spectrum  {X1#X2, . . . ,X  },  a correspond 
ing  eigenvector  and  generalised  eigenvector  set 

~XL 

h«l, 2, . . . ,m^ ; j«l, 2, . . . , k^ ; i«l, 2 , . . . , t } , and  a corresponding 
reciprocal  eigenvector  and  reciprocal  generalised  eigenvector 


set  {^C,b)  :c-l,2,...,mba;b-l,2,...,ka;a-l,2,...,t},  then 

A 

- ? («»> 

'iChW&Sl'ill'”  ’ i!^'3’  <27b) 


1. 


Li 


(m.  -l,j) 

‘^t-xi5+St9tSt>Jfx13  - ix,1 

(j"l»2, . • . , k^ ; i*l, 2 , . . . ,t) 


and 


♦x!'bl ' 


” Aa 


(27mjl) 


(28a) 
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where  is  the  hth  vector  in  the  jth  string  of  length 

-A. 

* /-  u\ 

associated  with  the  eigenvalue  and  where  ' 1 is 

A 

the  cth  vector  in  the  bth  string  of  length  m^  associated 

with  the  eigenvalue  Xft.  The  vectors  ( j*l# 2 , . . . ,k^) 

are  the  k^  eigenvectors  associated  with  the  eigenvalue  X^r 

whilst  the  remaining  vectors  in  each  of  the  k^  strings  of 

vectors  satisfying  equations  (27)  are  generalised  eigenvectors 

associated  with  the  eigenvalue  X..  Similarly,  the  vectors 
<ny  ,b) 

(b»l,2, . . . ,k  ) are  the  k reciprocal  eigenvectors 

A cl  A 

A 

associated  with  the  eigenvalue  X , whilst  the  remaining 
vectors  in  each  of  the  k strings  of  vectors  satisfying 
equations  (28)  are  reciprocal  generalised  eigenvectors  associated 
with  the  eigenvalue  X . The  total  number  of  vectors  associated 
with  the  eigenvalue  X^  in  each  set  is  evidently 


mf 


Z m 

g»l 


gf 


and 


(f»l,2,...,t)  (29) 


t 

n + pq  + i a Z mf  . (30) 

f-1  1 


Equations  (27)  and  (28)  can  be  written  in  the  form 
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where 


' A1 


Sx 


(c,b) 

a 


- Si5ii 


(c,b) 


(33) 


(34) 


and 


^,(a,b)x  (h,j) 
~Xa  ~Xi 


5 U 

ai  bj  ch 


(35) 


(h“  l(2r  • • . j“l»2|  • • • , k^ ; 1st  1 , 2 , « • . ft) 

(c-1,2, . . . *mba;b-l,2, . . . ,ka? a-1, 2 , . . . , t) 

Conversely,  if  equations  (31),  (32),  and  (35)  are  satisfied 
by  a self-con jugate  set  {X^, *2# • • • » *tI  of  complex  numbers 
and  corresponding  self-conjugate  sets  {xjh'^:  h-1,2 , . . . ,m. . t 

- J 1 

j-1, 2,  — ,k± ; i-1, 2,  — , t}  and  U<c'b)  sc-1,2, ...  .m^jb-l, 2, ...  ,ka#- 

a«l, 2, ..., t)  of  linearly  independent  vectors,  then  equations 

(33)  and  (34)  are  satisfied  by  an  (m+fc)  x(n+pq+i)  matrix 

such  that  {X^,Xj, . . . ,Xfc)  is  the  eigenvalue  spectrum  of  the 

matrix  (A^+B^G^C^)  , ik-1,2, . . . ,1^;  j-1,2, . . . ,k1» 

i“l,2,...,t}  constitutes  a corresponding  eigenvector  and 

generalised  eigenvector  set,  and  { 4>^Cr :c«l,2, . . . »n»baJ 

a 

b«l,2, . . . ,ka;a»l, 2, . . . ,t)  constitutes  a corresponding  reciprocal 
eigenvector  and  reciprocal  generalised  eigenvector  set.  It 
accordingly  follows  from  equations  (33)  and  (34)  respectively 
that  the  real  output-feedback  matrix 


'l  “ fel  * S?2  * • “p+pq+fl  ti£l  ' ?t^2  * *•*  * -i^p+pq+t^ 


-1 


(36) 
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and  the  real  transposed  output- feedback  matrix 

* CSl  ' $2  ' ***  ' Sm+i^  Q?i$l  ' ?l^2  ' •**  ' ?i$m+^ 


-1 


(37) 


assign  the  self-conjugate  eigenvalue  spectrum  {X^,!^  . . . ,Xfc} , 

the  corresponding  eigenvector  and  generalised  eigenvector  set 

{ X ^ ^ :h“l, 2 , . . . , ^ j j*li 2; • • • , k^ ; i»l ,2, • • • , t } / and  the 

corresponding  reciprocal  eigenvector  and  reciprocal  generalised 

eigenvector  set  {<t>*c'b)  :c=l,2, . . . ,mba;b«l,2, . . . ,kfi;a-l,2, . . . , t} 

a 

to  the  closed-loop  plant  matrix  (A^+B^G^C^)  in  case 
{ Xi , X2  > • • • 'Sj^p+pq+fc*  is  a subset  of  (p+pq+i)  linearly 
independent  members  of  the  set  (C.x|h' ^ ;h»l ,2 , . . . ,m. . ; 

••  A ^ J 1 

3®  1 > 2 , • • • f k^ ; i^l , 2 , • . • , t } and  a 

subset  of  (m+i)  linearly  independent  members  of  the  set 
{B'^c'b^  :c«l,2, . . . ,11^  ;b«l,2,  . . . ,ka;a=l,2,  . . . ,t}. 

d 

In  view  of  equations  (31),  (32),  (36),  and  (37),  the 
computation  of  is  thus  reduced  to  the  determination  (Porter 
and  D'Azzo  1978)  of  the  kernels  and  generalised  kernels  of 
each  of  the  t matrices 


- !*] 


(i»l,2, . . . ,t) 


(38) 


together  with  the  kernels  and  generalised  kernels  of  each  of 
the  t matrices 


( a*l ,2, •••ft) 


(39) 


followed  by  the  selection  of  sets  of  linearly  independent  self- 
conjugate vectors  (x{h'j)  sh-1,2, . . . ,mjl?  j-1,2, . . . fk^i-l,  2 , . . . , t} 
and  <*x  ' sc"^» • • • *n\j3a;b*l, 2 , . . . ,ka;a*l, 2, . . . ,t}  from  sub- 

A 

spaces  determined  by  the  kernels  and  generalised  kernels  of 
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i . 


§«,  d“l» 2'  * • • • t)  !i^a)  (aal/2,...,t),  respectively, 
such  that  the  orthogonality  conditions  (35)  are  satisfied. 

It  is  finally  evident  from  equations  (11),  (12),  and  (23) 
that  the  matrices  in  the  respective  state  and  output  equations 
of  the  required  ith-order  dynamic  compensator  are  determined 
by  the  sub-matrices  of  the  output-feedback  matrix  G. . 

J V 


jl 
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4.  Illustrative  Example 


The  results  presented  in  Sections  2 and  3 can  be  con- 
veniently illustrated  by  designing  an  error-actuated  dynamic 
controller  which  will  cause  the  output  of  the  controllable 
and  observable  linear  plant  governed  by  the  respective  state 
and  output  equations 


x^k+1) 
_x2  (k+1) 


0 , 1 


1 , -1/2J  [x2(k) 


x,  (k)1  fo]  fl 

+ u (k)  + d (k) 
x-(k)  1 -1 


y (k)  = |1  , 6] 


xx  (k) 
x2  ^ 


to  track  any  constant  command  input  r(k)  in  the  presence  of 
any  unmeasurable  constant  disturbance  input  d(k)  in  such  a 
way  that  e(k)  = 0 after  a finite  number  of  discrete-time 
intervals.  In  this  case  it  is  clear  that  r=s=q=l,  so 
that  the  open-loop  tracking  system  is  governed  by  the  respective 
state  and  output  equations 


x^ (k+1)  To  1 , Ol  f"x^ (k)~j  To' 

x9 (k+1)  = 1 , -1/2  ,0  x- (k)  + 1 u (k) 


z (k+1) 


-1  , 0 , 11  I z (k) 


‘o' 

'1* 

+ 

0 

r (k)  + 

-1 

d (k) 

(42) 

_1 

_0 

and 
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Y(k) 

z(k) 


1 * O 
0 , 0 


0* 

1 


x^k) 

x2(k) 

z(k)  J 


(43) 


This  system  is  clearly  controllable  and  observable  with 
vc  « 3 and  vQ  **  2 so  that  (Brasch  and  Pearson  1970)  equation 
(13)  indicates  that  l - 1.  Furthermore,  in  the  notation  of 
equations  (24)  , (25) , and  (26)  , it  follows  from  equations 
(42)  and  (43)  that 


(44) 


(45) 


(46) 


It  is  thus  evident  from  equations  (38) , (39) , (44)  , (45)  , 


and  (46)  that 
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?1(X) 


-A  i 0 , 0 ,0,0,0 

1 , -A-l/2  , 0 ,0,1,0 

-1  , 0 , 1-A  , 0 , 0 , 0 

_0  , 0 , 0 , “A  ,0,1 


and 


(47) 


T£  ( A) 


-A 

1 

0 

0 


f 


I 


/ 


t 


1 , -1  , o , 1 , 0 , 0~ 

-A-l/2  , 0 , O , O , 0 , 0 

0 , 1-A  , O , 0 , 1 , 0 

0 , 0 , —A  ,0,0,  1 


(48) 


In  order  to  design  an  error-actuated  dynamic  compensator  for 
the  open-loop  tracking  system  governed  by  equations  (42)  and 
(43)  such  that  the  error  is  eliminated  after  a finite  number 
of  discrete-time  intervals,  it  is  necessary  to  assign  the 
value  zero  to  all  the  eigenvalues  of  the  closed-loop  plant 
matrix.  Such  an  assignment  clearly  requires  that  t = 1 , 
m^  * 4,  and  A^  = 0 and  therefore,  in  consonance  with  the  results 
of  Rosenbrock  and  Hayton  (1977),  that  k^=  1 and  m1]L  = 4.  It 
follows  (Porter  and  D'Azzo  1978)  from  equation  (47)  that 


ker  S^O) 


span 


(49) 


and  similarly  (Porter  and  D'Azzo  1978)  from  equation  (48)  that 
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ker  T^(0) 


span 


1/2 

1 

0 

0 

-1 

0 

0 


0 

0 

1 

0 

1 

-1 

0 


(50) 


It  is  evident  from  equation  (49)  that  the  closed-loop  eigen- 
vector corresponding  to  the  eigenvalue  X1  = 0 must  be 
assigned  to  the  subspace 


Ex(0) 


span 


(51) 


whilst  the  associated  string  of  generalised  eigenvectors 

and  X^'1^  must  be  generated  in  accordance  with 

equations  (31)  , and  it  is  similarly  evident  from  equation  (50) 

(4  1) 

that  the  closed-loop  reciprocal  eigenvector  4>^  ' ' correspond- 
ing to  the  eigenvalue  = 0 must  be  assigned  to  the  subspace 


1^(0)  = span  j 


'1/2" 

o' 

"o' 

1 

0 

0 

0 

9 

1 

9 

0 

> 

_ 0 _ 

0_ 

JL_ 

(52) 


whilst  the  associated  string  of  generalised  reciprocal  eigen- 
vectors and  must  be  generated  in 


~o 


- 


1 I m H 


L8  - 


(32) . Since  the  vectors 


(53) 


(54) 


(55) 


(56) 
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are  clearly  such  that 

$oC,1)'Xoh,1)  " 6ch  (c,h«l,2,3,4)  (61) 

it  follows  from  equations  (36)  and  (37)  that  the  required 
eigenstructure  assignment  is  achieved  by  the  output  feed- 
back matrix 


-7/4  , 3/4  , 1/4' 
. 1/2  , -1/2  , -1/2 


(62) 


In  view  of  equations  (11) , (12) , (23) , and  (62)  , the  correspond- 
ing dynamic  compensator  for  the  open-loop  tracking  system 
governed  by  equations  (42)  and  (43)  is  governed  by  the 
respective  state  and  output  equations 

w(k+l)  = -1/2  w (k)  - 1/2  e (k)  - 1/2  z (k)  (63) 

and 

u(k)  ■ 1/3  w (k)  + 7/4  e (k)  + 3/4  z(k)  (64) 

so  that  the  required  error-actuated  dynamic  controller  is 
characterised  by  the  transfer  function 

T(z)  = u(z)/e(z)  = (7z2-z-2)/(4z+2)  (z-1)  . (65) 


It  can  be  readily  verified  that  the  closed-loop  tracking 
system  governed  by  equations  (42),  (43),  (63),  and  (64) 
tracks  any  constant  command  input  r(k)  in  the  presence  of 
any  constant  unmeasurable  disturbance  input  d(k)  in  such  a 
way  that 
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e (k)  = r (k)  - y(k)  = 0 (k»4,5,...)  . (66) 

4.  Conclusion 

In  this  paper,  the  method  of  entire  eigenstructure 
assignment  has  been  applied  to  the  design  of  linear  multi- 
variable  discrete-time  tracking  systems  incorporating  error- 
actuated  dynamic  controllers.  The  theory  developed  in  this 
paper  is  analogous  to  that  developed  by  Porter  and  Bradshaw 
(1978)  for  continuous-time  tracking  systems.  However,  in 
this  paper  the  theory  has  been  extended  in  order  to  allow  the 
assignment  of  confluent  eigenvalues  to  the  plant  matrix  of 
the  closed-loop  tracking  system.  It  is  therefore  possible  to 
apply  the  theory  to  the  design  of  error-actuated  dynamic 
controllers  which  eliminate  completely  the  error  between  the 
command  input  vector  and  the  output  vector  after  a finite 
number  of  discrete-time  intervals. 
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Abstract 

In  this  paper,  the  controllability  results  of  Bradshaw 
and  Porter  (1975)  are  applied  to  the  design  of  linear  mutli- 
variable  discrete-time  tracking  systems  incorporating  plants 
with  inaccessible  states  for  which  appropriate  error-actuated 
controllers  exist.  The  design  method  is  illustrated  by  the 
presentation  of  the  results  of  simulation  studies. 
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1.  Introduction 

It  Is  the  purpose  of  this  paper  to  develop  a basis 
for  the  design  of  linear  multivariable  discrete-time  tracking 
systems  incorporating  error-actuated  controllers  which  is 
directly  analogous  to  that  developed  for  continuous- time 
tracking  systems  by  Porter  and  Bradshaw  (1976) . The  control- 
lability results  of  Bradshaw  and  Porter  (1975)  are  accordingly 
applied  to  the  design  of  linear  multivariable  discrete- 
time tracking  systems  incorporating  plants  with  inaccessible 
states  for  which  appropriate  error-actuated  controllers  exist 
and  for  which  it  is  therefore  possible  to  achieve  satisfactory 
tracking  behaviour  without  the  need  to  incorporate  observers 
in  the  manner  of  Bradshaw  and  Porter  (1976) . Such  discrete- 
time tracking  systems  consist  of  a controllable  nth-order 
linear  multivariable  plant  governed  by  state  and  output 
equations  of  the  respective  forms 

x (k+1)  « Ax(k)  + Bu(k)  (1) 

and 

y (k)  - Cx(k)  (2) 

together  with  a controller  which  is  required  to  cause  the 
pxl  output  vector  y(k)  of  the  plant  to  track  a pxl  command 
input  vector  v(k)  in  the  sense  that 

lim  e(k)  » lim{v(k)  - y(k)>  * 0 (3) 

k-*»  ~ k-*<»> 

for  polynomial  command  inputs,  i.e.,  for  command  inputs  with 
the  property  that 
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A(r) v(k)  - 0 


where 


A(1)v(k) 
A(2)  v(k) 


v(k+l)  - v(k) 
A(1)v(k+1)  - A(1)v(k) 


A(m)v(k)  - A ^m”1)  v(k+l)  - A(m_1)v(k) 
*•«**«*» 


(4) 


It  is  important  to  note  that,  although  these  discrete-time 
tracking  systems  reduce  to  the  error-actuated  sampled-data 
servomechanisms  of  classical  control  theory  (Bergen  and 
Ragazzini  1954)  in  the  special  case  p ■ 1,  the  design  of 
error-actuated  multivariable  servomechanisms  in  the  general 
case  p > 1 is  always  non-trivial  - and  sometimes  impossible  - 
in  view  of  the  fact  that  the  assignment  of  prescribed  eigen- 
value spectra  by  error-actuated  controllers  and  by  output- 
feedback  controllers  (Kimura  1975)  are  essentially  equivalent. 

2.  Theory 

The  first  stage  in  the  design  of  the  required  error- 
actuated  controller  for  the  plant  governed  by  equations  (1) 
and  (2)  involves  the  introduction  (Bradshaw  and  Porter  1975) 
of  a vector  comparator  and  a series  of  r discrete-time  vector 
integrators  in  order  to  generate  the  r vectors  defined  by 
the  equations 
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z^k+1)  - z^k)  + e (k) 

z2(k+l)  - j2  (k)  + z^k) 
z3(k+l)  ■ z3(k)  + z2(k) 


zr(k+l)  - zr (k)  + zr_i (k) 


It  is  then  evident  from  equations  (1),  (2),  and  (6)  that  the 
open- loop  tracking  system  is  governed  by  a state  equation 
of  the  form 
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The  second  stage  in  the  design  of  the  error-actuated  controller 
involves  the  introduction  of  (r+1)  vector  feedback  loops  in 
order  to  generate  the  nutl  input  vector  u(k)  according  to  the 
error  control- law  equation 

I 

r 

u(k)  - K e (k)  + I K . z . (k)  , (8) 

~°~  i-l'1-1 

where  the  K. (i«0, 1,2, . . . ,r)  are  mxp  feedback  matrices.  It 
is  then  evident  from  equations  (7)  and  (8)  that  the  closed- 
loop  tracking  system  is  governed  by  a state  equation  of  the 
form 


( f 
| 


x(k+l) 

^(k+l) 

z2(k+l) 

z3(k+l) 


?r-l(k+1) 

sr(k+l) 


x(k) 

!l(k) 
z2  (k) 

z3  (k) 


5r(k) 


0 

«• 

0 

•• 


(9) 


a 


j 
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1 


u 

0 


In  the  case  of  pxl  vector  polynomial  command  inputs  of  the 
form 


v(k) 


la,  . k 
i-l~1_i 


i-1 


(10) 


it  is  clear  by  differencing  equation  (8)  (r-1)  times  that 
the  closed-loop  tracking  system  will  behave  so  that 

lim  A(r)z  (k)  - 0 (11) 

k-*» 


and  therefore,  in  view  of  equations  (6),  so  that  equation 
(3)  will  be  satisfied  if  the  error  control  law  (8)  can  be 
synthesised  in  such  a way  that  all  the  eigenvalues  of  the 
plant  matrix  of  the  closed- loop  tracking  system  governed  by 
equation  (9)  are  assigned  to  any  desired  locations  within 
the  unit  circle. 

However,  in  view  of  the  presence  of  the  sub-matrix 
(A-BKqC)  in  the  plant  matrix  of  the  closed-loop  tracking 
system,  an  error  control  law  of  this  class  will  not  always 
exist  (Kimura  1975)  even  if  the  open-loop  tracking  system 
governed  by  equation  (7)  is  controllable  in  the  sense  that 
(Bradshaw  and  Porter  1975) 


rank 


B 

0 


A-I 

-C 


n + p 


(12) 


It  is  nevertheless  evident  that  such  a control  law  will 
certainly  exist  if,  for  example,  a stabilising  state-feedback 
control  law  of  the  form 


I 

t 


j 
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r 

u(k)  = Kx(k)  + l K.z  (k)  (13) 

isl'1-1 

can  be  synthesised  such  that  there  exists  a matrix  K with 

-o 

the  special  property  that 

-KqC  = K . (14) 

However,  the  existence  of  a stabilising  error  control  law 
of  the  form  (8)  can  in  general  only  be  investigated 
systematically  by  using  decision  methods  in  the  manner  of 
Anderson,  Bose,  and  Jury  (1975) . 

3.  Illustrative  Example 


The  theory  presented  in  Section  2 can  be  conveniently 
illustrated  by  designing  an  error-actuated  controller  which 
will  cause  the  output  of  the  controllable  second-order  linear 
plant  governed  by  the  respective  state  and  output  equations 
(Bradshaw  and  Porter  1975) 


x^k+lf 

"0  , 1* 

~x1  (k)' 

ri  , in 

’u,  (kf 

3 

+ 

X 

_x2(k+l)_ 

-6  ' 5_ 

-x2  (k). 

i 

n 

O. 

-u2(k)- 

and 


Vi  00“ 

s 

'l  , o' 

'x1  (k)  " 

x2{k)~ 

-1  , 1_ 

*2  (k) 

to  track  the  command  input  vector 


i-* 

£ 
i 

3 

_2k 

v2(k)_ 

k 

(16) 


(0£k<«) 


(17) 
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In  this  case  it  is  clear  that  the  command  input  is  of  the 
form  (10)  with  r - 2,  so  that  the  state  equation  (7)  of  the 
open-loop  tracking  system  assumes  the  form 

0 , 1 , 0,0,  0,6 


“x^k+1) 

r 

x2 (k+1) 

Z11  (k+D 

z21(k+l) 

z12 (k+1) 

,z22(k+l). 

L 

1 , -1  , 0 , 1 , 0 , 0 

0 f 0 


’ 

’x^k) 

x2  (k) 

zlx (k) 

z21(M 

z12 

-z22(k)- 

i , r 
0 , 2 
0 , 0 
0 , 0 
0 , 0 
0 , o 


ux(k) 
•U2  <k> 


0 , 0 

0 , 0 

1 , 0 

o , 1 

0 , 0 

0 , 0 


2 k 
k 


rank 


1 , 1 , -1  , 1 

0 , 2 , -6  , 4 

0 , 0 , -1  , 0 

IP  , 0 , 1,-1 


- 4 


(18) 


where  z^tk)  is  the  ith  element  of  the  vector  z^ (k)  . Since 
(Bradshaw  and  Porter  1975) 


(19) 


the  controllability  condition  (12)  is  satisfied  in  this 
cases  it  is  therefore  certainly  possible  to  synthesise  a 
state-feedback  control  law  of  the  form  (13)  and  a correspond- 
ing error  control  law  of  the  form  (8)  such  that  the  eigenvalues 
of  the  plant  matrix  of  the  resulting  closed- loop  trackixxg 
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system  assume  arbitrary  values  since  also  the  output  matrix 
in  equation  (14)  is  invertible.  In  the  particular  case 
when  these  eigenvalues  are  all  assigned  the  value  zero  by 
the  implementation  of  the  error  control  law 


u^k) 

_u2 (k) 


2.5  , -2.5  , 1.5 
0.5  , 3.5  , 1.5 


-1.5  , 0.5 
1.5  , 0.5 


-0.5 

0.5 


e^k) 
e2  (k) 
*n(k) 

z2i(k) 

z12 (k) 
-z22 


(20) 


the  behaviour  of  the  initially  quiescent  tracking  system  is 
as  shown  by  the  full  lines  in  Figs  1 and  2 : it  is  evident 
from  Fig  1 that 


lim  e. (k)  - lim{v. (k)  - y. (k)  } » lirntv. (k)  - x.  (k)  } ■ 0 

A k-*-00  1 1 k-M»  1 


(21a) 


and  from  Fig  2 that 

lim  e- (k)  « lim{v, (k)  - y,(k)}  - lim{v0(k)  + x. (k)  - x« (k) } 
k—  2 k«-«  2 2 k>«  2 1 2 

(21b) 


as  required. 

The  corresponding  behaviour  of  the  initially  quiescent 
tracking  rystem  in  case  a state-feedback  control  law  is 
implemented  (Bradshaw  and  Porter  1975)  is  as  shown  by  the 
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dotted  lines  in  Figs  1 and  2:  it  is  again  evident  from  Figs  1 
and  2 that  equations  (21a)  and  (21b)  are  satisfied,  but  that 
the  transient  behaviour  of  the  tracking  system  incorporating 
a state-feedback  controller  (Bradshaw  and  Porter  1975)  is 
slower  and  less  oscillatory  than  the  corresponding  behaviour 
of  the  tracking  system  incorporating  an  error-actuated  controller. 

4.  Conclusions 

In  this  paper,  the  simple  matricial  methods  developed 
by  Bradshaw  and  Porter  (1975)  for  the  design  of  linear 
multivariable  discrete- time  tracking  systems  for  plants  with 
accessible  states  have  been  applied  to  the  design  of  linear 
multivariable  discrete-time  tracking  systems  incorporating 
plants  with  inaccessible  states  for  which  appropriate  error- 
actuated  controllers  exist.  The  results  of  simulation  studies 
have  been  presented  which  indicate  that  the  transient  behaviour 
of  tracking  systems  of  the  latter  class  is  faster  but  more 
oscillatory  than  the  corresponding  behaviour  of  tracking 
systems  incorporating  state-feedback  controllers. 
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Li 


Abstract 


It  is  known  (Porter  and  Bradshaw  1978a)  that,  in  the 
case  of  self-conjugate  distinct  eigenvalue  spectra,  the 
closed-loop  eigenstructure  assignable  by  output  feedback  is 
constrained  by  the  requirement  that  the  eigenvectors  and 
reciprocal  eigenvectors  lie  in  well-defined  subspaces.  In 
this  paper,  a technique  is  presented  which  can  be  used  to 
select  the  eigenvectors  and  reciprocal  eigenvectors  from 
these  subspaces  in  the  case  of  appropriately  augmented 
(Kimura  1975)  controllable  and  observable  continuous -time 
systems.  This  technique  is  ideally  suited  to  digital- 
computer  implementation  and  therefore  greatly  facilitates 
the  synthesis  of  both  static  (Porter  and  Bradshaw  1978a) 
and  dynamic  (Porter  and  Bradshaw  1978b)  output- feedback 


controllers. 


1.  Introduction 


It  has  bean  shown  (Porter  and  Bradshaw  1978a, b)  that  the 
method  of  entire  eigenstructure  assignment  can  be  applied  to 
the  design  of  output-feedback  controllers  for  multivariable 
linear  continuous -time  systems  governed  by  state  and  output 
equations  of  the  respective  forms 

x(t)  ■ Ax(t)  + Bu(t)  (1) 

**»  «*»  «•  •%» 

and 

^(t)  - Cx  ( t)  , (2) 

where  A G ^ . bER115™,  c£  Rpxn,  rank  B = m,  and  rank 

**•  «N«  «« 

C ■ p.  Thus,  if  output  feedback  Is  applied  to  the  system 
governed  by  the  state  equation  (1)  in  accordance  with  the 
control- law  equation 

u ( t ) « Gy (t)  (3) 

««•  *v 

and  the  output- feedback  matrix  G £ SUCh  that  the 

closed-loop  plant  matrix  (A+BGC)  has  the  self-conjugate 

distinct  eigenvalue  spectrum  { X. , X0 , . . . , X >,  then  the  corre- 

i i n 

sponding  eigenvector  and  reciprocal  eigenvector  sets  (u, , u- , . . . ,u  } 

- 1 -2  - n 

and  {v. , v, , . . . , v ) must  be  such  that 
-i  -i  ~n 


Sfj 


E kerfA-Xjl  , B]  (j-1,2, . ..  ,n) 


(4) 


?i 


E ker[A'-XiI  , C']  (i-1,2, . . . ,n)  , (5) 


J 


y 


and 
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ri“j  ■ {ij 


(i»  j“li 2» • • ■ «n)  • (6) 


The  output- feedback  matrix  is  then  given  by  the  formulae 


9 " fel  ' *2  ' * * * * i 


» CUj  » • • • » Cul 


I"1 


5'  ■ [*i  > *2 • s'y2 ?'yj 


where,  in  this  paper,  the  state  and  output  equations  (1)  and 

(2)  represent  appropriately  augmented  (Kimura  1975)  con- 
trollable and  observable  continuous-time  systems. 

Thus,  the  synthesis  of  the  output- feedback  control  law 

(3)  requires  the  selection  of  linearly  independent  sets  of 
vectors  {u^Uj, . . . ,un)  an<*  ^v^»v2' • * • » vn^  ^rom  the  respec- 
tive subspaces  defined  by  relations  (4)  and  (5)  which  satisfy 
the  orthogonality  condition  (6) . It  is  shown  that  this 
selection  can  be  effected  by  performing  restricted  elementary 
row  and  column  operations  on  matrices  formed  from  the  spanning 
vectors  of  these  subspaces,  and  that  the  resulting  synthesis 
procedure  is  therefore  ideally  suited  to  digital  computer 
implementation . 


2,  Theory 


The  first  stage  in  the  synthesis  of  the  control  law  is 
clearly  the  determination  (Porter  and  D’Azzo  1977)  of  the 
n kernels 
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S(Xj)  ■ ker(_A-XjI  , B]  « span 


IVvl 

J :k-l, . . . ,m  • 

HwJ  J 


( j"l*  2 , . , . , n) 


and  tha  n kernels 


'•(»1)  “ ker(A'-XI  , C']  ■ span 


!k(V 

£k(V 


• k* 1 f • • • t p 


( i®  1 f 2f  • » • ^ n) 


It  follows  from  relation  (4)  and  equation  (9)  that 


UjEU(Xj)  « span  {Xk(Xj)!k-l,...,m) 


( j»l,2, . . . ,n) 


and  from  relation  (5)  and  equation  (10)  that 


YiE  - span  {^(X^  ik-1 p} 


(i-l,2,...,n)  , (12) 


that  is 


u,  - X(X  )n. 
■* J “ J -J 


( j“l » 2 , . . . ,n) 


5i!'<xi> 


(i-l,2,...,n)  , (14) 


where 


?(Xj)  " ^l(Xj)  ' Xj  ( X j ) 


^m(  Xj)1 


( j"l»2, . . . ,n) 


(15) 
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and 

4(X^)  m » $2^1^  ' ***  ' -p^i^  * 

(i«l,2,...,n)  (16) 

where  n^  is  an  mxl  vector  and  i'  is  an  lxp  vector.  Equations 
(13)  and  (14)  can  be  more  conveniently  expressed  in  the  form 

Uj  - XU^JNj  (j«l,2,. ,.,n)  (17) 

and 


■ L^4*(X^)  ( iml , 2 , . . • , n)  , (18) 

where  u.  and  n.  are  the  first  columns  of  the  nxm  matrix  U. 
and  the  mxm  matrix  N ^ respectively,  and  v^  and  are  the 
first  rows  of  the  pxn  matrix  and  the  pxp  matrix 
respectively.  The  orthogonality  condition  (7)  then  requires 
that  the  element  in  the  first  row  and  the  first  column 
of  each  of  the  pxm  matrices 

5ij  - 


(i,j-l,2,...,n)  (19) 

be  such  that 

m^j  ■ ^ij  ( i t j"l» 2 , • * . , n)  , (20) 

where 

° ^ *4#(X^)X(Xj)  (i»  j*l»  2 , . . . ,n)  • (21) 

It  is  evident  that  the  condition  (20)  can  be  satisfied  by 
performing  restricted  elementary  row  and  column  operations 


n 
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on  the  npxnm  matrix 
M«0>  - 

according  to  the  equation 


(i,3«l,2,...,n)  (22) 


M - LM(o)N 


(23) 


where  the  npxnm  matrix 


*?  - C?13] 


the  npxnp  matrix 


( 1 » j-1, 2 , . . . ,n) , (24) 


L - diag[L.^] 


( 1—1 , 2 / . • • ,n)  , (25) 


and  the  nmxnm  matrix 


N - diag[Nj] 


(3*1,2, . • • ,n) 


(26) 


These  computations  can  be  conveniently  organized  in 
the  following  steps  in  view  of  the  results  of  Kimura  (1975)  : 


(1)  Set  M(o)  » [M^J]  , L - Inp  , and  N 


(o) 


^nm' 


(2)  By  restricted  elementary  column  operations  on  and 

Nj°^  (1*1,2, . . . ,n$  j-1, 2, . . . ,n-m)  determine  u^£[  U(X^) 

( j"l» 2, . . . ,n-m)  such  that  {u, ,u-,...,u  } is  a linearly 

-i  -a  -n-m  * 

independent  set,  -►  N*1*  , and 

*■***»  *• 


M(1>  - 


(27) 


(3)  By  restricted  elementary  row  operations  on  m| ^ and 
(i*n-m+l,n-m+2, . . . ,n; j*l,2, . . . ,n)  determine 

(i»n-m+l,n-m+2, . . . ,n)  such  that  vju^  ■ m|2) 
- 0 (i*n-m+l,n-m+2, . . . ,n; j*l,2, . . . ,n-m) , M(1)  m(2)  , 


(28) 


M 


(2) 


L(1)M<°)N(1) 


(2) 


(4)  By  restricted  elementary  column  operations  on 


(1) 


and  (i»l,2, . . . ,n; j«n-m+l,n-m+2 , . . . ,n)  determine 

UjE  ( j»n-m+l,n-m+2 , . . . ,n)  such  that  v'u^  ■ mj3^ 


* 0 (i=n-m+l,n-m+2 , . . . ,n; j*n-m+l,n-m+2, . . . ,n;i^j) , 
M(2)  ■*  M(3)  , N(1) 


N(2)  , and 


m<3)  = l(1)m(o)n(2) 


(29) 


,(3) 


(5)  By  restricted  elementary  row  operations  on  ^ and 
L f ^ (i=l, 2 , . . . ,n-m;  j*l,2,...,n)  determine  v|£i/'(^) 


(4) 

(i=l, 2, . . . ,n-m)  such  that  vf u . = m . . = 0 (i=l,2,...,r 

.1-j  Xj 


j*n-m+l,n-m+2 , . . . ,n) , M 


(3) 


M<4)  , LU) 


i/  2 ^ , and 


M 


(4) 


l(2)m(o)n(2) 


(30) 


(6)  By  restricted  elementary  row  or  column  operations  on 


M 


(4) 


(i=l,2, ..  . ,n; j=l#  2, ...  ,n)  normalize  v'  or 


such  that  v'ui  - mj3)  * 1 (i-1 ,2, . . . ,n) , M(4)  - M(5) 


(2) 


L(3)  , N(2) 


N(3)  , and 


M 


(5) 


L(3)M(o)N(3) 


M 


LM(o)N 


(31) 


(7)  Compute  the  u^  ( j-1,2, . . . ,n)  using  equation  (17)  and 
compute  the  v(  (i«l,2, . . . ,n)  using  equation  (18). 
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In  certain  pathological  cases  (Kimura  1975) , special 
spectra  exist  for  which  no  corresponding  output  feedback 

matrix  exists  and  for  which  this  computational  procedure 

. 

therefore  fails:  in  such  cases , however,  it  is  only  necessary 
slightly  to  perturb  the  spectra  in  order  to  obtain  solutions. 
It  is  also  possible  for  this  computational  procedure  to 
fail  for  certain  pathological  choices  of  u.£li(X.) 

( j = l, 2, . . . ,n-m)  in  step  (2):  in  such  cases,  however,  it 
is  only  necessary  slightly  to  perturb  the  u^  £ ti  ( X ^ ) 

( j=l,2, . . . ,n-m) . 


3.  Illustrative  Example  l 


The  procedure  can  be  conveniently  illustrated  by  the 
synthesis  of  an  output- feedback  control  law  for  the  continuous- 
time system  governed  by  the  respective  state  and  output 
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U(-l) 


1 


span 


1 


9 


U(-2) 


span 


4 


-1 


0 

0 


U(-3) 


and 


U(-4) 


0 

1 

* ) 

’ 

-3 

o' 

9 

0 

span  ■ 

-1 

9 

0 

w 

0 

1 

• 

-4' 

o' 

16 

0 

span 

-1 

9 

0 

0 

1 

9 


9 


(43) 


(44) 


(45) 


(46) 


and  it  is  similarly  evident  from  equations  (12),  (39),  (40), 
(41),  and  (42)  that  the  closed-loop  reciprocal  eigenvectors 
corresponding  to  the  eigenvalue  spectrum  (34)  must  be  assigned 
to  the  respective  subspaces 


l/(-l)  ■ span 


-2' 

0 

* 

1 

0 

0 

0 

9 

1 

9 

0 

_0_ 

0 

1 

(47) 
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Then,  by  performing  the  restricted  elementary  column  operations 
corresponding  to  setting 


?1 


(52) 


and 


n 


2 


(53) 


it  follows  from  equation  (27)  that 


(54) 


from  equation  (28)  that 
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diag 


1/3, 1/3, 1/3} 

-9/2,-ll,-5|  4, 

13  ,5  J-7/6, -14/3, -5/3’ 

0 , 1 , 0 { 

1 

•h 

o 

o 

H 

o 

1 ,0  S 0 , 1 ,0 

j 

. o ,-14,  1 ! 

0 -14,  1 ! 0, 

-7/5,1  ! 0 , -7/5  , 1 

(59) 


and 


N * diag 

'5  , 0 J 1 #0 

1 

1 

0 J 

1 

1 / o' 

-7,lj-7,l 

-14 

, 1 i 

-23  , 1. 

. (60) 


It  therefore  follows  from  equation  (17)  that 


fr 


{Yl'Y2'Y3'Y4} 


-5' 

-2" 

' -3' 

' -4’ 

5 

4 

9 

16 

-5 

9 

-1 

9 

-1 

9 

-1 

-7 

-7 

-14 

-23 

(61) 


and  from  equation  (18)  that 


<Yl'Y2'Y3'Y4} 


“2/3' 

'27/2' 

-16' 

' 35/6  ' 

1/3 

-9/2 

4 

-7/6 

1/3 

9 

-11 

9 

13 

9 

-14/3 

1/3. 

-5 

5 

(62) 


In  view  of  equations  (4)  and  (9)  the  results  (61)  imply 
that 


-5' 

-10' 

'-33' 

-76' 

9 

9 

9 

7 

. 14 

42 

l J 

92 

L J 

{Yl'Y2'Y3'Y4}  " 1 I - I ' I I ' I I » I l)»  <63> 

and  in  view  of  equations  (5)  and  (10)  the  results  (62)  imply 


that 
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{?1'?2'53'?4} 


-1/3  , 

-1/3 


-67/2 


-161/6' 
56/3 
. 20/3 


It  finally  follows  from  either  equation  (7)  or  equation  (8) 
that  the  output- feedback  control  law  (Porter  and  Bradshaw 
1978b) 


u ( t) 


-47  , 34  , 10 

49  , -35  , -11 


assigns  the  spectrum  (34),  the  eigenvectors  (61),  and  the 
reciprocal  eigenvectors  (62)  to  the  closed-loop  plant  matrix 
of  the  system  governed  by  equations  (32)  and  (33). 


4.  Conclusion 


It  is  known  (Porter  and  Bradshaw  1978a)  that,  in  the 
case  of  self-conjugate  distinct  eigenvalue  spectra,  the 
closed-loop  eigenstructure  assignable  by  output  feedback  is 
constrained  by  the  requirement  that  the  eigenvectors  and 
reciprocal  eigenvectors  lie  in  well-defined  subspaces.  In 
this  paper,  a technique  has  been  presented  which  can  be  used 
to  select  the  eigenvectors  and  reciprocal  eigenvectors  from 
these  subspaces  in  the  case  of  appropriately  augmented 
(Kimura  1975)  controllable  and  observable  continuous-time 
systems  by  performing  restricted  elementary  row  and  column 
operations  on  matrices  formed  from  the  spanning  vectors  of 
these  subspaces.  This  technique  is  ideally  suited  to  digital- 
computer  implementation  and  therefore  greatly  facilitates  the 


154 


synthesis  of  both  static  (Porter  and  Bradshaw  1978a)  and 
dynamic  (Porter  and  Bradshaw  1978b)  output- feedback  controllers. 
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Abstract 

In  view  of  the  fundamental  new  insights  into  the 
. structure  of  linear  multivariable  continuous-time  systems 
provided  by  the  method  of  entire  eigen structure  assignment, 
the  design  of  dynamic  compensators  is  equivalent  to  the 
selection  of  pairwise-orthogonal  eigenvectors  and  reciprocal 
eigenvectors  from  two  families  of  well-defined  subspaces 
which  are  parametrised  by  associated  self-conjugate  eigenvalue 
spectra.  This  selection  is  effected  by  the  use  of  a powerful 
new  algorithm  which  requires  the  performance  of  restricted 
elementary  row  and  column  operations  on  matrices  formed  from 
the  spanning  vectors  of  these  subspaces.  The  digital  computer 
implementation  of  the  resulting  procedure  incorporating  this 
algorithm  is  described  and  is  illustrated  by  the  design  of  an 
error-actuated  dynamic  compensator  for  a linear  multivariable 
plant. 


!' 


1,  Introduction 


In  most  practical  cases,  it  is  of  course  impossible  co 
implement  state-feedback  control  laws  since  the  state  of  the 
plant  is  inaccessible  and  only  the  plant  output  is  available 
for  control  purposes.  The  method  of  entire  eigenstructure 
assignment  (Porter  and  D'Azzo,  1977)  has  accordingly  been 
applied  to  the  design  of  linear  multivariable  continuous- 
time output- feedback  regulators  by  Porter  and  Bradshaw  (1978a) 
It  has  been  shown  that,  in  the  case  of  self-conjugate  distinct 
eigenvalue  spectra,  the  closed-loop  eigenstructure  assign- 
able by  output  feedback  is  constrained  by  the  requirement 
that  the  eigenvectors  and  reciprocal  eigenvectors  of  the 
closed-loop  plant  matrix  lie  in  two  families  of  well-defined 
subspaces  and  satisfy  appropriate  orthogonality  conditions. 

In  contrast,  the  closed-loop  eigenstructure  assignable  by 
state  feedback  (Moore,  1976)  is  constrained  only  by  the 
requirement  that  the  eigenvectors  of  the  closed- loop  plant 
matrix  lie  in  just  one  family  of  well-defined  subspaces.  It 
is  because  of  the  severe  constraints  on  the  closed-loop 
eigenstructure  assignable  by  output  feedback  that  it  is 
frequently  impossible  to  achieve  satisfactory  closed-loop 
behaviour  by  means  of  static  continuous-time  output-feedback 
regulators,  and  that  it  is  consequently  necessary  to  introduce 
dynamic  compensators  (Brash  and  Pearson,  1970;  Kimura,  1975) . 
However,  it  has  been  shown  by  Porter  and  Bradshaw  (1978b) 
that  the  design  of  such  dynamic  compensators  can  be  effected 
by  applying  the  method  of  entire  eigenstructure  assignment  to 
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appropriately  augmented  (Kimura,  1975)  continuous-time 
systems . 

In  view  of  these  fundamental  new  insights  into  the 
structure  of  linear  multivariable  systems,  the  design  of 
dynamic  compensators  is  equivalent  to  the  selection  of 
pairwise-orthogonal  eigenvectors  and  reciprocal  eigenvectors 
from  two  families  of  well-defined  subspaces  which  are 
parametrised  by  associated  self-conjugate  eigenvalue  spectra. 
This  selection  can  be  effected  by  the  use  of  a powerful  new 
algorithm  (Bradshaw,  Fletcher,  and  Porter,  1978)  which 
requires  the  performance  of  restricted  elementary  row  and 
column  operations  on  matrices  formed  from  the  spanning  vectors 
of  these  subspaces.  The  digital  computer  implementation  of 
a procedure  incorporating  this  algorithm  is  described  and 
is  illustrated  by  the  design  of  an  error-actuated  dynamic 
compensator  for  a linear  multivariable  plant.  The  principal 
computational  attraction  of  the  procedure  is  that  no  operations 
with  polynomial  matrices  are  involved,  so  that  error-actuated 
dynamic  compensators  for  large-scale  systems  can  be  readily 
designed. 

2.  Compensator  Structure 

The  linear  multivariable  continuous-time  tracking 
systems  considered  by  Porter  and  Bradshaw  (1978b)  consist 
of  a controllable  and  observable  nth-order  plant  governed 
by  state  and  output  equations  of  the  respective  forms 


x(t)  = Ax  ( t)  + Bu  ( t)  + Dd  ( t) 


(1) 
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y (t)  - Cx(t)  , (2) 

where  x(t)  E Rn,  u(t)  E Y (t)  E RP,  d(t)  E ^ A E RnXI\ 
B ERnXm»  C ERPXn»  D ERnXh*  rank  B - m,  and  rank  C - p, 
together  with  an  error-actuated  dynamic  compensator  which 
is  required  to  cause  the  output  vector,  y(t) , to  track  a 
command  input  vector,  r(t),  in  the  sense  that 


lim  e(t)  * lim{r(t) -y (t) } = 0 

£•*■«  f*100 


for  unmeasurable  command  and  disturbance  inputs  of  the  respective 
forms 


E ct.  , t 
i-1  1 1 


d(t)  - e e ,tJ 
i-i  1 1 


Such  an  error-actuated  dynamic  compensator  is  governed  by 
state  and  output  equations  of  the  respective  forms 


w ( t)  - Fw ( t)  + Ge ( t)  + E H.z.(t) 

i-1  1 1 


Kw  ( t)  + Le  ( t)  + E M.  z.  (t) 

i-1  1 1 


where 


q - max(r,s)  , (9) 


w(t)ER*»  e(t)E*P*  z±(t)ERP  (i-l,2,...,q)  , u(t)  E Rm/ 

f£RW(  g£R1xP,  H1£l!lxp  (i«l,2,...,q)  , 

L E Rmxp , (i=l»2,  . . . ,q)  , and  (Kimura,  1975) 

l =■  max(0,n-m-p+l)  . (10) 


It  is  then  evident  from  equations  (1) , (2) , (6) , (7) , and 
(8)  that  the  closed- loop  system  is  governed  by  state  and 
output  equations  of  the  respective  forms 
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’ y (t)  ' 

C#  0 i • • • i 0 t 0 

zL(t) 

0 , I , « • • , 0 , 0 

r 

SB 

ZqU) 

0,  0 , • • • , ip  , o 

w(t) 

p , 0 , ...  , 0 , I, 

x(t) 

Z^t) 


2q(t) 

w(t) 


(12) 


In  view  of  equations  (4)  and  (5),  it  is  clear  by  differentiating 
equation  (11)  (q-1)  times  that  the  closed-loop  system  will 

behave  so  that 


lim 

t-*«0 


0 


(13) 


and  therefore,  in  view  of  equations  (8) , so  that  equation 
(3)  will  be  satisfied  if  the  error-actuated  dynamic 
compensator  governed  by  equations  (6)  and  (7)  is  designed 
such  that  all  the  eigenvalues  of  the  plant  matrix  of  the 
closed-loop  system  governed  by  equations  (11)  and  (12)  are 
assigned  to  the  open  left-half  of  the  complex  plane. 

It  is  evident  from  equations  (11)  and  (12)  that  such 
a compensator  can  be  designed  by  the  synthesis  of  an  appro- 
priate output-feedback  control  law  of  the  form 

Ui(t)  “ (14) 

for  the  augmented  open- loop  system  governed  by  state  and 
output  equations  of  the  respective  forms 

xt(t)  - A^U)  + (15) 


and 
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B , 0 
0 , 0 


and 


(22) 


(23) 


Thus,  if  the  (m+l)  x(p+pq+£)  output- feedback  matrix  G^  is 
such  that  the  closed- loop  plant  matrix  (A^+B^G^C^)  has  the 
self-conjugate  distinct  eigenvalue  spectrum  A = {X^,X2,..., 
Xn+pq+^},  then  the  corresponding  eigenvector  and  reciprocal 

eigenvector  sets  (u^Uj un+pq+Jl>  and  (v^Vj 

must  be  such  that 


Uj 

-Wj- 


G 


ker[A1-Xj1n+pq+ll 


B4J  (j-1,2 


,n+pq+i)  (24) 


ker[Aj-AiIn+pq+jl 


q]  (i-1,2 


,n+pq+fc)  (25) 


and 


viuj  * su 


(i, j-1,2, . . . ,n+pq+6) . (26) 


The  output- feedback  matrix  is  then  given  by  the  equivalent 
formulas 


Gfi.  " £wi  ' W2  ' * * * ' Wp+pq+*l  ^C£U1  ' C£U2  ' * ' * ' CaUp+pq+J 


-1 


(27) 


and 

Q[  - [*!  » *2  ' • 


# *m+J  ^BJvi  ' B)tv2  ' •**  * B£vm+J 


-1 


(28) 


In  view  of  equations  (24),  (25),  (27),  and  (28),  the 
computation  of  G£  is  reduced  to  the  determination  of  the 
kernels  of  each  of  the  n+pq+fc  matrices 

Si(V  “ ^Ai"XjI  ' BJ  (j-1,2 n+pq+4)  (29) 

together  with  the  kernels  of  each  of  the  n+pq+s,  matrices 

Tf(Xi)  * fAi"XiI  ' CJ1  (i-l,2,...,n+pq+f)  (30) 

followed  by  the  selection  of  sets  of  linearly  independent 
self-conjugate  vectors  {u1#u2, . . . ,uR  +jt)  and  {v^Vj,..., 
vn+pg+£>  from  subspaces  determined  by  the  kernels  of  S^(X^) 

( j"l,2, , . . ,n+pq+f)  and  T^(X^)  (i«*l,2  , . . . ,n+pq+£)  , respectively, 
such  that  the  orthogonality  conditions  (26)  are  satisfied. 

It  is  finally  evident  from  equations  (6) , (7) , and  (20)  that 
the  matrices  in  the  respective  state  and  output  equations 
of  the  required  ith-order  error-actuated  dynamic  compensator 
are  determined  by  the  sub-matrices  of  the  output-feedback 


V 


matrix  G 
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3.  Compensator  Design  Procedure 

3.1  System  Augmentation  Procedure 

The  first  stage  in  the  compensator  design  procedure  involves 
the  formation  of  the  augmented  plant,  input,  and  output 
matrices  by  the  following  steps  which  constitute  the  routine 
AUGMENT : 

(i)  Set  q = max(r,s) ; 

(ii)  Set  ! ® max (0,n-m-p+l) 

(iii)  Form  the  augmented  open-loop  plant,  input,  and 
output  matrices  A^,  B^,  and  C^. 

3.2  Kernel  Computation  Procedure 

The  second  stage  in  the  compensator  design  procedure  involves 
the  computation  of  the  closed-loop  eigenvector  and  reciprocal 
eigenvector  subspaces  by  the  following  steps  which  constitute 
the  routine  KERNELS: 


(i)  Select  the  closed- loop  eigenvalue  spectrum  A = 
{A1,A2,...,An+pq+1}; 


(ii) 

Form  S^Uj) 

- (vy 

( j»l,2, . . . ,n+pq+i) ; 

(iii) 

Form  T'U^ 

- 

(i"l,2, . . . ,n+pq+l)  ; 

(iv) 

Compute 

ker  S & ( Aj ) 


span 


IVV" 

k‘V- 


:k=l , 2 


,m+l  ) 


( j“l,2 , . . . ,n+pq+fc) ; 
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(v)  Compute 


ker  TJUJ 


span 


w 

w 


:k-l,2 


. ,p+pq+i  « 


(i«l,2, . . . ,n+pq+i) ; 

(vi)  Form  X(X^)  - » X2<Xj>  ' • xm+i*Xj^ 

( j“l»2, . . . ,n+pq+i)  i 

( vii )  Form  Q(Xj)  “ r ^2  (Xj)  / • • • » ( x j ) ] 

( j-lf  2, . . . ,n+pq+t)  ; 

(viii)  Form  *(X1)  - [4>1(Xi)  , ^(X^  f ...  , ♦p+pq+t(xi)] 

(i»l,2# . . . ,n+pq+l) ; 

(ix)  Form  Z(X^)  *■  [;1(xi)  > 52(xj_)  ' ***  ' ^p+pq+t^i^ 

(i»l,2 , . . . ,n+pq+i)  . 


3.3  Eigenvector  Selection  Procedure 

The  third  stage  in  the  compensator  design  procedure  involves 
the  selection  of  the  pairwise-orthogonal  closed-loop 
eigenvectors  and  reciprocal  eigenvectors  from  the  respective 
subspaces  im  X(X^)  ( j-1,2 , . . . ,n+pq+£)  and  im  MX^) 

(i«l,2 , . . . ,n+pq+fc)  by  the  following  steps  which  constitute 
the  routine  SELECT: 

\L)  Select  u^  £ im  X(X^)  ( j»l,2  , . . n+pq-m)  such  that 
{u^,Uj » . • • »un-fpq_m}  is  a linearly  independent  set? 


1 


i 


1 
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(ii)  Compute  vi  £ im  4>  ( X^)  (i=n+pq-m+l,n+pq-m+2  , . . . ,n+pq+£) 
such  that  vfu^  = 0 (i=n+pq-m+l  ,n+pq-m+2  , . . . /n+pq+i.; 
j=»l » 2 , . . . ,n+pq-m)  ; 

(iii)  Compute  u^  £im  X ( X ^ ) ( j=n+pq-m+l,n+pq-m+2  , . . . ,n+pq+fl) 
such  that  v^Uj  = 0 (i=l,2 , . . . ,n+pq-m; j=n+pq-m+l, 
n+pq-m+2 , . . . ,n+pq+4.) ; 

(iv)  Compute  v^^  £ im  ^ ( X ^ ) (i=l , 2 , . . . ,n+pq-m)  such  that 
v^Uj  = O (i=l,2, . . . ,n+pq-m; j=n+pq-m+l,n+pq-m+2, . . . , 
n+pq+t) ; 

(v)  Normalise  vi  or  such  that  v'^  = 1 (i=l , 2 , . . . ,n+pq-Ht) 
3.4  Compensator  Matrix  Computation  Procedure 


The  final  stage  in  the  compensator  design  procedure  involves 
the  computation  of  the  compensator  matrices  by  the  following 
steps  which  constitute  the  routine  COMPENSATE: 


(i)  Select  a set  {C^  , C£u2  , ...  , Vp+pq+J,}  of 

linearly  independent  vectors  and  a set  {B£v1  , b^v2 
B'vm+^}  of  linearly  independent  vectors; 


(ii)  Compute  the  output- feedback  matrix 


GJt  “ ^Wl,W2,,**'Wp+pq+J  ^CJ.Ul,CilU2,***'C)lUp+pq+^ 


-1 


and  the  transposed  output-feedback  matrix 


Gl  " [>i'z2'---'znHj  tB£vl'B£v2'*  * ‘ 


-1 


(iii)  Form  the  compensator  matrices  K#LfM^,M2 


F Hj / • • • f Hg • 


i • • • t cUid 
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4.  Illustrative  Example 


This  procedure  can  be  conveniently  illustrated  by 
designing  an  error-actuated  dynamic  compensator  which  will 
cause  the  output  of  the  controllable  and  observable  linear 
plant  governed  by  the  respective  state  and  output  equations 


and 


x1  ( t) 

r 

x2  (t) 

n 

x3  ( t) 

x4(t)_ 

0 , 1,0,0 

0 , -1  , -1  , 0 

1 , 0 , 0 , 0 


Xj.tt)  ' 

x2  (t) 

x3  (t) 

“ ■ » 

4J 

X 

1 

H* 

o 

' 0 ' 

o 

o 

u^t) 

+ 

0 

o 

o 

-u2(t)- 

1 

1 

o 

H* 

-1 

d(t) 


(31) 


*<l 

H* 

ft 

i 

1 , -i  , 2 , O' 

.y2(t). 

p , 0 ,1,1 

xx(t) 

x2  (t) 
x3(t) 
Lx4(t) 


(32) 


to  track  any  constant  command  input  r(t)  = [r^t)  , r2(t)]'  £ R4 
in  the  presence  of  any  unmeasurable  constant  disturbance  input 
d(t)  R. 


In  this  case,  the  outputs  of  the  routines  AUGMENT, 
KERNELS,  SELECT,  and  COMPENSATE  when  A = {-1.0, -1 . 5 ,-2 .0 , 
-2. 5, -3. 0,-3. 5, -4.0}  are  listed  in  the  Appendix.  These 
listings  indicate  that  the  required  error-actuated  dynamic 
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0 


compensator  is  governed  by  the  respective  state  and  output 
equations 


w(t)  - -1.989  w ( t)  + [-0.5431  , -2.258] 


+ [-1.260  , 5.227] 


ex(t) 
1*2  (t> 


z^t) 

.*2(t) 


(33) 


and 


u.  (tf 

'-19.40  ‘ 

'10.26  , -51.98  ' 

X 

= 

w(t)  + 

1 

.u2 

1.823 

. 0.6437  , 4.961 

— i 

CM 

0) 

I 

-9.521  , 0.2316 

2.645  , -3.509 


zx(t) 

z2(t). 


(34) 


where  [e^t)  , e2(t)]'  - [r^tj-y^t)  , r2 ( t)  -y2 (t)]  ' £ R2 
and  [z  x ( t ) , Zj  ( t)  ] ■ [e^t)  , e2(t)]'^R2. 

5,  Conclusion 

The  method  of  entire  eigenstructure  assignment  has 
yielded  fundamental  new  insights  into  the  structure  of  linear 
multivariable  systems  and,  in  particular,  into  the  closed- 
loop  eigenstructure  assignable  by  output  feedback  (Porter 
and  Bradshaw,  1978a, b) . The  design  of  dynamic  compensators 
has  accordingly  been  reduced  to  the  selection  (Bradshaw, 
Fletcher,  and  Porter,  1978)  of  pairwise-orthogonal  eigenvectors 
and  reciprocal  eigenvectors  from  two  families  of  well-defined 
subspaces  which  are  parametrised  by  associated  self-conjugate 
eigenvalue  spectra.  The  resulting  procedure  for  the  design 


of  dynamic  compensators  is  computationally  attractive  since 
its  constituent  routines  AUGMENT , KERNELS , SELECT , and 
COMPENSATE  involve  only  numerically  stable  operations. 

Indeed,  the  entire  procedure  has  been  coded  in  FORTRAN  for 
the  routine  computer-aided  design  of  error-actuated  dynamic 
compensators,  and  forms  part  of  a comprehensive  suite  of 
design  procedures  for  various  classes  of  controllers  for  both 
continuous -time  and  discrete-time  linear  multivariable 
systems . 


Acknowledgment 

This  research  was  supported  in  part  by  the  Wright 

Aeronautical  Laboratories,  United  States  Air  Force  under 

Grant  AF0SR-76-3005B. 

References 

Bradshaw,  A,  L R Fletcher,  and  B Porter  (1978) . Synthesis 
of  output-feedback  control  laws  for  linear  multivariable 
continuous -time  systems.  Int  J Syst  Sci,  9^  1331-1340. 

Brasch,  F M and  J B Pearson  (1970).  Pole-placement  using 
dynamic  compensators.  IEEE  Trans  Autom  Control,  AC-15, 
34-43. 

Kimura,  H (1975).  Pole  assignment  by  gain  output  feedback. 
IEEE  Trans  Autom  Control,  AC-20,  509-516. 


y 


171  - 


Moore,  B C (1976).  On  the  flexibility  offered  by  state 
feedback  in  multivariable  systems  beyond  closed-loop 
eigenvalue  assignment.  IEEE  Trans  Autom  Control,  AC-21 , 
689-692. 

Porter,  B and  A Bradshaw  (1978) . Design  of  linear  multi- 
variable  continuous -time  output-feedback  regulators.  Int 
J Syst  Sci,  445-450. 

Porter,  B and  A Bradshaw  (1978).  Design  of , linear  multi- 
variable  continuous-time  tracking  systems  incorporating 
error-actuated  dynamic  controllers.  Int  J Syst  Sci,  9, 
627-637. 

Porter,  B and  J J D'Azzo  (1977) . Algorithm  for  the  synthesis 
of  state-feedback  regulators  by  entire  eigenstructure 
assignment.  Electron  Lett,  13,  230-231. 


E i! 

• 

172  - 

I U ■ 

Appendix 

1 L 

AU0M6NI 

AUGMENTED  MATRIX  A 

0 OOOOC 

OO 

I.  OOOOC 

00 

0.  OOOOC 

oo 

0 OOOOC 

oo 

0 OOOOC 

OO 

0. 

OOOOE 

OO 

O.  OOOOE 

00 

0 OOOOC 

oo 

-1.  OOOOC 

oo 

-1.  OOOOC 

oo 

O OOOOC 

OO 

0.  OOOOE 

00 

0. 

OOOOE 

OO 

0 OOOOC 

oo 

I OOOOC 

OO 

0.  OOOOC 

00 

0.  OOOOC 

oo 

0.  OOOOE 

oo 

o ooooe 

00 

o. 

ooooe 

oo 

O.  OOOOC 

00 

-1  OOOOC 

oo 

0.  OOOOC 

00 

I.  OOOOE 

oo 

2.  OOOOE 

00 

0.  OOOOC 

oo 

0. 

OOOOE 

00 

o ooooe 

00 

-I  OOOOC 

00 

1 OOOOC 

oo 

-2.  OOOOC 

oo 

O.  OOOOE 

00 

0.  OOOOE 

00 

0. 

OOOOE 

00 

O.  OOOOC 

00 

0 OOOOC 

00 

O.  OOOOC 

oo 

OOOOC 

OO 

-1.  OOOOE 

oo 

0 OOOOE 

00 

0. 

OOOOC 

oo 

0.  OOOOC 

00 

0 OOOOC 

oo 

0.  OOOOC 

oo 

a ooooe 

oo 

O.  OOOOE 

oo 

o.  ooooe 

oo 

0. 

OOOOE 

00 

0.  OOOOE 

00 

AUGMENTED  MATRIX  B 

* 

1 OOOOC 

00 

0.  OOOOC 

oo 

o.  ooooe 

00 

0 OOOOC 

00 

0 OOOOC 

00 

0.  OOOOC 

oo 

0 OOOOC 

oo 

I.  OOOOE 

oo 

a ooooe 

00 

0 OOOOC 

oo 

0.  OOOOC 

00 

O.  OOOOC 

oo 

0 OOOOC 

oo 

0 OOOOC 

oo 

0.  OOOOC 

00 

0 OOOOC 

00 

0 OOOOC 

00 

0.  OOOOE 

00 

0.  OOOOC 

00 

0.  OOOOC 

oo 

I.  OOOOC 

oo 

OIUMCNTCII  MATRIX  C 

1.  OOOOC 

oo 

-1.  OOOOC 

00 

2.  OOOOC 

00 

0.  OOOOE 

oo 

0.  OOOOE 

oo 

0 

OOOOE 

00 

0 OOOOE 

00 

0.  OOOOC 

oo 

O.  OOOOC 

oo 

I.  ooooe 

00 

1.  OOOOC 

oo 

0.  OOOOE 

00 

0. 

ooooe 

00 

0 OOOOC 

00 

0.  OOOOC 

00 

0.  OOOOC 

oo 

o.  ooooe 

00 

1 OOOOC 

oo 

0 OOOOE 

00 

0 

OOOOE 

00 

0 OOOOE 

00 

0.  OOOOC 

00 

I.  OOOOC 

oo 

0.  OOOOC 

o 

0.  OOOOC 

oo 

o ooooe 

00 

0 

OOOOE 

00 

0 OOOOE 

oo 

0 OOOOC 

oo 

0.  OOOOC 

oo 

I.  OOOOE 

CO 

0.  OOOOC 

00 

0 OOOOC 

oo 

0. 

OOOOE 

00 

0.  OOOOE 

00 

1 i 

KERNE!  S 

UFPt-H  SPANNING  VECTORS 

5 CHKLAMHUA) 

1 OOOOC 

oo 

0 OOOOC 

oo 

0.  OOOOC 

00 

1.  OOOOE 

oo 

0.  OOOOC 

oo 

0. 

ooooe 

00 

1.  OOOOE 

00 

0 OOOOC 

00 

1.  OOOOE 

00 

0.  OOOOC 

00 

0 OOOOE 

00 

1.  OOOOC 

oo 

0. 

OOOOC 

oo 

0.  OOOOE 

00 

0 OOOOE 

00 

0 OOOOC 

00 

0.  OOOOE 

00 

0 OOOOE 

oo 

5.  ooooe-oi 

0. 

OOOOE 

00 

0 OOOOE 

00 

0.  OOOOE 

00 

0 OOOOC 

00 

0 OOOOC 

oo 

2 8370E-01 

-1.  429OE-01 

0. 

OOOOC 

00 

2 3OOOE-01 

S I OOOOE 

oo 

-I  OOOOE 

00 

o ooooe 

00 

6.  6670E-01 

-3  <.380C-12 

0. 

ooooe 

oo 

3.  OOOOE-OI 

0 OOOOE 

00 

0 OOOOE 

00 

0.  OOOOE 

00 

1 9030E-01 

2 3810E-0I 

0. 

OOOOC 

00 

1.  2SOOE-01 

0.  OOOOE 

00 

0 OOOOE 

oo 

1.  OOOOC 

oo 

0.  OOOOC 

oo 

0.  OOOOE 

oo 

1. 

OOOOC 

00 

0.  OOOOE 

00 

0 OOOOE 

00 

0 OOOOE 

oo 

1 OOOOE 

00 

0 OOOOE 

00 

0.  OOOOE 

00 

I. 

OOOOE 

oo 

0 OOOOE 

oo 

1.  OOOOE 

00 

0 OOOOE 

00 

0 OOOOE 

00 

1 OOOOE 

00 

0 OOOOE 

00 

0. 

OOOOE 

oo 

1.  OOOOE 

00 

1.  00O0F 

00 

0 OOOOC 

00 

-3  C380E- I 2 

1 OOOOE 

00 

0 OOOOE 

00 

0 

OOOOE 

00 

2.  OOOOE 

00 

-2  SOOOE- 

■ot 

0 OOOOE 

00 

2 2;'20C-0! 

0 OOOOE 

00 

0 OOOOE 

00 

2 

OOOOE-OI 

-4.  OOOOE-OI 

\ 3 0000E-01 

0 OOOOE 

00 

4 OOOOE- 

-01 

0.  OOOOE 

oc 

0 OOOOE 

oo 

3. 

3330E-0I 

I.  OOOOE 

00 

| 3 7300E-01 

0.  OOOOE 

00 

8 88*0E-02 

0.  OOOOE 

00 

0.  OOOOE 

00 

6. 

6A70C-02 

3.  3330E- 

-01 

0.  OOOOE 

oo 

I OOOOE 

00 

0 OOOOE 

00 

0 OOOOC 

00 

1.  OOOOE 

00 

0 

OOOOE 

oo 

0.  OOOOE 

oo 

0 OOOOE 

00 

1 OOOOC 

00 

0.  OOOOE 

00 

0 OOOOE 

00 

I.  OOOOC 

00 

0. 

OOOOE 

00 

0.  OOOOE 

00 

0 OOOOE 

oo 

0 OOOOE 

00 

1.  OOOOE 

00 

0.  OOOOE 

00 

0 OOOOE 

00 

1. 

OOOOE 

00 

0 OOOOC 

00 

0 OOOOE 

00 

0 OOOOE 

00 

2.  SOOOE 

00 

0 OOOOC 

00 

3.  6380E-12 

3. 

OOOOE 

00 

0.  OOOOC 

oo 

0.  OOOOE 

00 

1 8180E-01 

-4  S450E-01 

o ooooe 

00 

1.  6&70E-01 

-3.  OOOOE -Ol 

0 OOOOE 

00 

f 0.  OOOOE 

00 

2 8S70E- 

-ot 

t.  1430E 

00 

0.  OOOOE 

00 

2.  3000E-01 

1. 

2SOOE 

00 

0 OOOOC 

oo 

0 OOOOE 

oo 

3 1930E- 

-02 

3.  8440E-01 

0 OOOOE 

oo 

4 1670E- 

-02 

6. 

2SOOE-0 1 

o ooooe 

00 

1.  OOOOE 

00 

0 OOOOE 

00 

0.  OOOOC 

00 

1.  OOOOE 

00 

0 OOOOC 

oo 

0. 

OOOOE 

00 

1.  OOOOE 

oo 

I 1 OWER  SPANNING  VICTORS 

i OMCOA(LAMHHA) 

£ - I OOOOE 

oo 

-1.  OOOOC 

oo 

0 OOOOC 

00 

-1.  SOOOE 

00 

-1  OOOOE 

00 

0. 

OOOOE 

00 

-2  OOOOE 

00 

I -1.  OOOOE 

00 

0.  OOOOE 

00 

0 OOOOE 

00 

-1.  OOOOE 

00 

-7  SOOOE -01 

0 

OOOOE 

oo 

-l  OOOOC 

oo 

i 0 OOOOE 

00 

0 OOOOE 

00 

-1  OOOOE 

00 

0.  OOOOE 

00 

0 OOOOC 

00 

-1. 

SOOOE 

oo 

0.  OOOOE 

oo 

I |*  -1  OOOOE 

00 

0 OOOOC 

00 

-2  SOOOE 

oo 

-1.  OOOOE 

00 

0 OOOOE 

00 

-3. 

OOOOE 

00 

-i.  ooooe 

00 

S -2  OOOOE 

oo 

0 OOOOC 

oo 

-1  OOOOE 

00 

-3  7300C 

00 

0 OOOOE 

00 

-1 

OOOOE 

oo 

-C.  OOOOE 

oo 

p L 0.  OOOOE 

00 

-2  OOOOC 

00 

0 OOOOC 

00 

0 OOOOE 

00 

-2  SOOOE 

oo 

0 

OOOOC 

oo 

0.  OOOOC 

00 

| 0 OOOOE 

00 

- 3 3OO0E 

00 

-1.  OOOOC 

00 

0 OOOOE 

oo 

-4  OOOOE 

00 

-I. 

OOOOC 

oo 

0 OOOOC 

oo 

H 0 OOOOE 

oo 

-I  OOOOC 

oo 

-e  7oooc 

00 

0 OOOOE 

00 

-1  OOOOC 

00 

-1 

2000C 

ot 

o ooooe 

00 

l\\  -3.  OOOOC 

oo 

0 OOOOC 

00 

0.  OOOOC 

00 

-3.  30006 

00 

0 OOOOC 

00 

0. 

00  X* 

oo 

0.  OOOOE 

00 

173 


I*  HER  REC  1 I KOCAL  SPANNING  VECTORS 
PHI (LAMMOA) 


1. 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

00 

0 OOOOE 

00 

0 

OOOOE 

00 

1 

OOOOE 

00 

0 OOOOE 

oo 

0. 

ooooe 

00 

1 OOOOE 

00 

0. 

OOOOE 

00 

0 OOOOE 

00 

0. 

ooooe 

00 

0 

ooooe 

00 

1.  OOOOE 

oo-. 

-2. 

ooooe 

00 

-1. 

OOOOE 

00 

0 

OOOOE 

00 

0 OOOOE 

00 

0. 

OOOOE 

00 

-4. 

2300E 

oo 

-1.  2SOOE 

<sr 

0 

ooooe 

00 

-1 

OOOOE 

oo 

0 

OOOOE 

00 

0.  OOOOE 

00 

0 

OOOOE 

00 

-1. 

7300E 

00 

-7.  3000C-01 

0. 

OOOOE 

00 

0. 

OOOOE 

oo 

t. 

OOOOE 

00 

0 OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

oo 

0 OOOOE 

00 

0 

OOOOE 

oo 

0. 

OOOOE 

00 

0. 

OOOOC 

oo 

1 OOOOE 

oo 

0. 

ooooe 

00 

0. 

OOOOE 

oo 

0.  OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

o 

OOOOC 

CO 

0.  OOOOE 

oo 

1. 

OOOOE 

00 

0. 

OOOOE 

00 

o ooooe 

00 

0. 

ooooc 

oo 

0. 

OOOOC 

oo 

0. 

ooooe 

00 

1.  OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

o ooooe 

oo 

0. 

ooooe 

oo 

0. 

OOOOE 

00 

0. 

ooooc 

00 

0.  OOOOE 

00 

ft. 

OOOOE 

00 

0. 

OOOOE 

00 

o ooooe 

00 

0 

IHKKJE 

00 

0. 

ooooc 

00 

0. 

ooooe 

00 

-7.  OOOOE 

oo 

-2 

ooooe 

oo 

0. 

OOOOE 

00 

0 OOOOE 

00 

o.  ooooe 

00 

0. 

OOOOE 

00 

0 

ooooc 

00 

-4  OOOOE 

00 

-1. 

ooooe 

00 

0. 

OOOOE 

00 

0 OOOOE 

00 

i 

OOOOE 

00 

0. 

ooooe 

00 

0 

ooooe 

00 

0.  OOOOE 

00 

0. 

OOOOE 

00 

1. 

OOOOE 

00 

0 OOOOE 

00 

0 

ooooc 

oo 

1. 

OOOOE 

00 

0 

OOOOE 

00 

0 OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

1 OOOOE 

oo 

0 

OOOOE 

oo 

0. 

OOOOE 

00 

1. 

OOOOE 

00 

o ooooe 

00 

0 

OOOOE 

oo 

0 

ooooe 

00 

0 OOOOE 

00 

0. 

OOOOE 

oo 

1. 

OOOOE 

oo 

0 

ooooe 

oo 

0.  OOOOE 

00 

0. 

ooooe 

00 

0. 

OOOOE 

00 

1 OOOOE 

oo 

0. 

OOOOE 

00 

0 

OOOOE 

00 

1.  ooooc 

00 

0 OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0 OOOOE 

oo 

0. 

OOOOC 

oo 

-1. 

02S0E 

01 

-3. 

2300C 

00 

0 OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

00 

-l  4000E 

01 

0 

OOOOE 

oo 

-4 

7SOOE 

00 

-1.  7300E 

00 

0.  OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

oo 

-1.  OOOOE 

01 

0 

OOOOC 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

1.  OOOOE 

00 

0. 

ooooe 

00 

0 

OOOOE 

00 

0.  OOOOE 

00 

0. 

ooooc 

00 

0. 

OOOOE 

oo 

0. 

OOOOE 

00 

0 OOOOE 

00 

1. 

OOOOE 

oo 

0. 

OOOOE 

oo 

0.  OOOOE 

00 

1. 

ooooc 

00 

0. 

OOOOE 

oo 

o. 

OOOOE 

oo 

0.  OOOOE 

oo 

0. 

OOOOE 

00 

1. 

OOOOE 

oo 

0.  OOOOE 

00 

0. 

ooooc 

oo 

0 

OOOOE 

00 

0. 

ooooe 

oo 

o ooooe 

oo 

1 

OOOOE 

00 

0 

OOOOE 

oo 

0.  OOOOE 

00 

1. 

ooooc 

oo 

0 

OOOOE 

oo 

0 

ooooc 

oo 

0.  OOOOE 

00 

o. 

OOOOE 

00 

1.  OOOOE 

00 

o ooooe 

00 

-3  OOOOC 

oo 

0. 

OOOOE 

00 

0. 

ooooc 

oo 

0.  OOOOE 

00 

-1.  B230E 

01 

-7. 

2300E 

00 

0.  OOOOE 

00 

-3.  OOOU* 

oo 

0. 

OOOOE 

00 

0 

ooooe 

oo 

0.  OOOOE 

oo 

-1. 

3730E 

01 

-4 

7300E 

00 

0.  OOOOE 

00 

0 

ooooc 

00 

t. 

OOOOE 

oo 

0 

OOOOE 

00 

0 OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

1.  OOOOE 

00 

0 

ooooc 

00 

0. 

OOOOE 

00 

1. 

OOOOE 

00 

0.  OOOOE 

00 

0. 

OOOOE 

oo 

0 

OOOoE 

00 

0 OOOOE 

00 

0 

ooooc 

oo 

0 

OOOOE 

00 

0 

ooooe 

oo 

i.  ooooe 

oo 

0 

ooooe 

00 

0. 

ooooe 

00 

0 OOOOE 

00 

0 

ooooc 

oo 

o. 

OOOOE 

00 

1. 

OOOOE 

00 

0.  OOOOE 

00 

0. 

ooooe 

00 

0. 

OOOOE 

00 

0.  OOOOE 

00 

0. 

ooooc 

00 

0. 

OOOOE 

oo 

0. 

OOOOE 

00 

1 OOOOE 

00 

0. 

ooooe 

00 

0. 

OOOOE 

00 

0 OOOOE 

00 

0 

ooooc 

00 

0. 

OOOOE 

00 

-2 

30O0E 

01 

-1  OOOOE 

01 

0. 

OOOOE 

00 

o. 

OOOOE 

00 

0 OOOOE 

00 

0 

ooooc 

oo 

0. 

OOOOE 

00 

-i  ooooe 

01 

-7  OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0 OOOOE 

00 

0. 

ooooe 

oo 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

0.  OOOOE 

00 

1 

ooooc 

oo 

0. 

OOOOE 

00 

0 OOOOE 

00 

1. 

OOOOE 

00 

0. 

OOOOE 

oo 

0. 

OOOOE 

00 

0 OOOOE 

oo 

0. 

ooooc 

oo 

1. 

OOOOE 

00 

0 OOOOE 

00 

0. 

ooooc 

oo 

I. 

OOOOE 

00 

0. 

OOOOE 

00 

0 OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

1 OOOOC 

oo 

I MWER  RECIPROCAL  SPANNING  VECTORS 
2ETAU  AtIBDA) 


t. 

ooooc 

OO 

-3. 

4370E-12 

1. 

OOOOE 

oo 

-3. 

6360E-12 

0. 

OOOOE 

OO 

1. 

ooooc 

OO 

5 OOOOC-Ol 

0 

ooooe 

00 

3. 

OOOOE 

00 

0 

OOOOE 

00 

1. 

OOOOE 

oo 

0. 

OOOOE 

oo 

6. 

1230C 

OO 

2. 

6230C 

00 

0. 

ooooc 

00 

0. 

OOOOE 

00 

-1. 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

00 

-3. 

3420C-II 

-e. 

8C80E-I2 

0 

ooooe 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

-1 

OOOOE 

00 

0 

OOOOE 

00 

0 

ooooc 

OO 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

00 

-1. 

OOOOE 

00 

0. 

ooooc 

00 

0 

OOOOE 

00 

1. 

OOOOE 

00 

-1. 

0910E-1 1 

0. 

OOOOE 

00 

i.  ooooe 

00 

1 

OOOOE 

00 

1. 

OOOOE 

00 

-1. 

0910E-1 1 

0. 

OOOOE 

00 

1 

ooooe 

00 

0. 

OOOOE 

00 

1. 

6000E 

01 

4 

OOOOE 

00 

0 

OOOOE 

00 

1 

OOOOE 

00 

-1 

oOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

r 

ooooe 

00 

0. 

OOOOE 

00 

-2 

OOOOE 

00 

0. 

OOOOE 

00 

0. 

ooooe 

00 

-1. 

3000E 

00 

0. 

OOOOE 

oo 

0 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

oo 

-2. 

OOOOE 

00 

0. 

OOOOE 

oo 

0 

OOOOE 

00 

-1. 

3000E 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0. 

ooooe 

00 

1. 

OOOOE 

00 

1. 

3000E 

00 

t 

ooooe 

1)0 

0. 

OOOOE 

00 

0 

OOOOE 

00 

1. 

OOOOE 

00 

0 

OOOOE 

00 

3. 

0370E 

01 

7. 

073OE 

oo 

0 

ooooe 

00 

1. 

OOOOE 

00 

0. 

OOOOE 

00 

3.  OOOOE 

01 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

-2 

3000E 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

oo 

0. 

OOOOE 

oo 

0. 

OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

-2. 

SOOOE 

00 

0. 

OOOOE 

00 

0. 

ooooe 

00 

-2. 

OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

-2 

30OOE 

oo 

0.  ooooc 

00 

2. 

OOOOE 

00 

1 

OOOOE 

oo 

0 

OOOOE 

00 

0. 

OOOOE 

00 

1 

OOOOE 

00 

2. 

aoooE 

oo 

1. 

ooooc 

oo 

1. 

S.OOOC 

01 

0 

OOOOE 

00 

1. 

OOOOE 

00 

0 

OOOOE 

00 

7 

3430E 

01 

2 

4130E 

01 

0. 

ooooc 

oo 

0. 

OOOOE 

00 

-3 

OOOOE 

00 

0 

OOOOE 

00 

0 

ooooe 

00 

0 

OOOOE 

00 

0. 

OOOOE 

oo 

-3. 

SOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

00 

-3. 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

oo 

0 

OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOE 

00 

-3 

OOOOE 

00 

b 

OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

-7 

V300E-13 

0 

OOOOE 

00 

1 

OOOOE 

00 

3 

OOOOE 

00 

t. 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

00 

I. 

OOOOE 

00 

0 

OOOOE 

00 

1 

obooe 

02 

4 

2000E 

01 

0 

OOOOE 

00 

1 

ooooe 

00 

0 

OOOOE 

00 

0. 

00001 

00 

0 

OOOoE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

-4 

OOOOE 

00 

0 

OOOOE 

00 

0. 

OOOOC 

00 

-3. 

soooe 

00 

0 

OOOOE 

00 

0. 

ooooc 

00 

0 

OOOOE 

00 

0. 

OOOOC 

00 

-4. 

OOOOE 

00 

0 

OOOOE 

00 

0. 

ooooc 

00 

-3 

3000E 

00 

0. 

OOOOE 

00 

0. 

OOOOE 

00 

0 

OOOOE 

00 

0 

OOOOE 

oo 

-4 

OOOOE 

00 

nm 


- 174  - 


SELECT 

CLOT'-tll-LOOf* 

ui.n 

t OOOOE  oo 
1 OOOOE  00 
0 00006  00 
3 33306-01 
0 00006  00 
3 3 <306-01 
I CK100E  00 


EIGENVECTORS 

t.  OOOOE  00 

1.  OOOOE  00 

5 OOOOE-Ot 
1.  42706-01 

6 66706-01 
4 8260E-01 
1 OOOOE  00 


t.  OOOOE  00 
0 00006  00 
0 OOOOE  00 

2.  50006-01 
5 00006-01 
1.  25006-01 
0 00006  00 


I.  OOOOE  00 
0 OOOOE  00 
-3  *3806-12 
2 2720E-01 

4.  OOOOE-Ol 
8.  88706-02 
0 OOOOE  00 


-1  07406  02 

1 OOOOE  00 

2 OOOOE  00 
-2  22706  01 
-3.  54506  01 
-6  75706  00 

2 41006"  OO 


-2.  7070E  01 

3.  38806-01 
8 4700E-01 
-5.  43706  00 
-7  71806  00 
-1.  3120E  OO 
1 OOOOE  00 


l.  OOOOE  00 
-1.  22306-02 
-3  6*806-02 

1 72806-01 

2 34706-01 

3 40206-02 
-4  12406-02 


Cl  USED-LOOT- 
V<  J> 

l 70606-01 
-2  12706  00 

1 76506  00 

2 13206  00 
-3  78306  00 

7.  45*06  00 
-2  50206-01 


HEClWlOCflL  EIGENVECTORS 

-4.  64506-01  -7  43506  OO  6 16306  01 

5.  27406  00  6.  5S70E  01  -4.  27206  02 

-4  64406  00  -7.  71106  01  7.  4040E  02 

-3.  15806  00  -3.  5840E  01  3 2780E  02 

6.  1770E  OO  6.  8160E  01  -7  00706  02 

-t  4680E  01  -1  3350E  02  1 6650E  03 

I 1I50E  00  -1.  6760E  00  -3.  0200E  02 


1 

25406 

00 

-3. 

70706 

00 

-2 

36106 

01 

-7. 

40406 

00 

2. 

08606 

01 

1. 

50106 

02 

1 7460E 

01 

-8. 

35306 

01 

-7 

57806 

02 

7 6730E 

OO 

-4. 

80706 

01 

-6. 

25606 

02 

-1. 

78206 

01 

7. 

21806 

01 

7. 

7710E 

02 

4. 

1710E 

01 

-1. 

62706 

02 

-1. 

67606 

03 

-1. 

1050E 

01 

5. 

3180E 

01 

6. 

40806 

02 

COMFENSO I E 

OUT  KIT-FEE  UBACK  MATRIX  0 

-»  02606  01  5 1V806  01  -7  52106  00  2 31606-01  -1.  7400E  01 

-6  43706-01  -4  7610E  00  2.  64506  00  -3  50706  00  1.  8230E  00 

5 41306-01  -2  25806  00  -1.  26006  00  5 2370E  00  -1.  7870E  00 


